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I INTOODLCTION 


A.  General 

When  an  aircraft  or  other  flight  vehicle  is  operated  in  precipita- 
tion containing  ice  crystals,  frictional  electrification  associated  with 

particle  impact  causes  the  impinging  particles  to  acquire  a net  charge 

1-5* 

and  to  deposit  an  equal  and  opposite  charge  on  the  vehicle.  As 

Shu vn  in  Figure  1,  charging  occurs  on  frontal  metallic  s'  ■ dielectric 

6,7 

portions  of  the  vehicle.  Although  the  charge  deposited  by  a single 


FIGURE  1 FRICTIONAL  CHARGING  OF  AIRCRAFT  BY  IWPINGING  PARTICLES 


crystal  changes  the  potential  of  the  aircraft  only  slightly  (of  the 

order  of  0,01  volt  change  lor  the  case  of  a KC-135  struck  by  a cirr is- 
l 

cloud  crystal  ) , the  particle  impact  rate  in  a typical  cloud  is  sufficient 

that  the  vehicle  reaches  a potential  of  hundreds  of  kilovolts  in  less 
•I 

than  a second. 

Charge  deposited  on  dielectric  surfaces  such  as  the  radome  or 
canopy  is  bound  there  because  these  surfaces  are  insulators.  As  a 

♦References  are  listed  •'it  the  end  of  the  report, 
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result,  under  precipitation  charging  conditions,  it  is  possible  for  a 
potential  difference  of  hundreds  of  kilovolts  to  exist  between  a dielec- 
tric surface  and  the  surrounding  airframe  structure. 

The  electrification  of  the  vehicle  is  of  relatively  little  conse- 
quence in  itself  because  the  energies  involved  are  small,  and  since  the 
electrostatic  fields  do  not  penetrate  to  the  interior.  It  is  the  conse- 
quences of  the  electrification  (Figure  2)  that  ar  ’ of  concern.  When  the 

CORONA 

DISCHARGES 

FROM 

EXTREMITIES 


FIGURE  2 NOISE  SOURCES  ASSOCIATED  WITH  STATIC  ELECTRIFICATION 

’'Chicle  potential  reaches  roughly  100  kV,  the  electric  field  intensity 

at  the  aircraft  extremities  becomes  sufficiently  high  that  electrical 

8 

breakdown  of  the  air  (corona  discharge)  occurs.  At  aircraft  operating 

altitudes,  the  corona  breakdown  from  the  extremities  occurs  not  as  a 

continuous  flow  of  charge  but  as  a series  of  pulses  with  roughly  10  ns 

I’ise  time  and  200  ns  duration,  and  it  therefore  generates  radio  noise 

■1 , , 8 

over  a broad  spectrum. 

A similar  situation  exists  on  the  dielectric  frontal  surfaces.  As 
charge  continues  to  accumulate  on  the  dielectric,  the  potential  to  the 
airframe  rises  until  the  electric-field  intensity  at  the  dielectric 
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surface  becomes  sufficiently  high  that  voltage  breakdown  (streamer  dis- 
charge) occurs  across  the  plastic  surface.  A surface  streamer  involves 

the  rapid  transfer  of  charge  over  a substantial  distance,  and  also 

6,7 

generates  serious  radio  frequency  interference. 

The  degree  to  which  the  radio  frequency  noise  generated  by  corona 

and  streamer  discharges  couples  into  electronic  systems  on  the  flight 

vehicle  is  determined  by  the  relative  locations  of  the  noise  source,  and 

the  "antenna"  via  which  the  noise  is  coupled  into  the  affected  system. 

In  addition,  the  coupling  depends  on  frequency,  the  size  of  the  vehicle, 

4 5 7 

and  the  size  of  the  antenna. 

B.  Scope  of  Present  Program 

On  earlier  programs,  various  aspects  of  the  problem  of  precipita- 
tion-static noise  generation  and  coupling  were  studied  analytically  and 
experimentally  both  in  the  laboratory  and  in  flight  tests.  Althou[n  the 
earlier  work  was  directed  largely  toward  specific  aircraft,  the  results 
are  capable  of  being  generalized  and  applied  to  similar  vehicles.  Un- 
fortunately, the  results  are  spread  over  a large  number  of  reports,  each 
of  which  treats  a limited  part  of  the  overall  problem.  Thus,  the  EMC 
engineer  is  in  the  position  of  having  to  be  familiar  with  all  of  the 
publications  in  considerable  depth  if  he  is  to  apply  the  results  of  the 
earlier  work  to  his  problems. 

The  purpose  of  the  present  program  is  to  systematize  the  results 
of  earlier  precipitation-static  studies,  and  to  prepare  mathematical 
models  based  on  this  earlier  work  to  permit  the  results  to  be  incor- 
porated in  the  Air  Force  Intrasystem  EMC  Analysis  Program  currently 
9 

under  development. 

It  was  recognized  at  the  outset  that  the  scope  of  the  program 
would  be  circiun .scribed  by  e.xisting  data  on  noise  sources  and  coupling. 
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It  was  noted,  however,  that  existing  results  could  be  applied  to  a con- 
siderable variety  of  situations,  and  that  incorporating  a precipitation- 
static  model  into  the  Intrasystem  EMC  Analysis  Program  at  this  time 
would  be  of  considerable  immediate  benefit  to  the  EMC  engineer.  Further 
refinements  of  the  precipitation-static  model  can  be  accomplished  as 
new  experimental  data  become  available  and  as  new  analytical  techniques 
are  developed. 

In  preparing  the  precipitation  static  model  (available  as  \'ol  II 
of  this  report),  one  of  the  objectives  was  that  it  be  compatible  with 
the  Air  Force  Intrasystem  EMC  Analysis  Program  ( lAP) . As  work  progressed, 
it  was  determined  that  the  structure  and  functioning  of  the  precipitation 
static  model  was  such  that  it  lent  Itself  to  embodiment  as  a "stand-alone" 
program.  Since  it  appeared  likely  that  there  would  be  further  refinemtuit 
and  development  of  the  lAP,  it  was  decided  that  the  precipitation  static 
model  should  indeed  be  generated  in  the  form  of  a stand-alone  program 
that  can  be  called  when  needed.  The  outputs  of  the  precipitation  static 
model  are  values  of  precipitation  static  noise  field  intensity  at  selected 
locations  on  the  aircraft.  These  form  natural  input  data  lor  the  lAP. 


4 


II  PRECIPITATION  CHARGING,  NOISE  GENERATION,  AND  COUPLING 


A.  General 

As  indicated  earlier,  considerable  work  has  been  done  on  the  in- 
vestigation of  precipitation-slatic  noise  generation  and  coupling 
processes.  The  essential  results  of  this  work  will  be  presented  here 
together  with  an  indication  of  the  special  considerations  or  limitations 
that  might  apply  to  the  result  used.  An  effort  will  be  made  to  avoid 
repeating  derivations  and  analyses  here.  Instead,  extensive  reference 
will  be  made  to  the  original  work„ 

L.  Noise  Sources 

1.  Corona  Discharge 

In  Ref.  8 it  was  shown  that  the  individual  pulses  associated 
with  a corona  discharge  can  be  well  approximated  by  a decaying  ex- 
ponential with  zero  rise  time,  as  follows: 


f (t^  = A e~^^  (1) 

whore  A is  the  pulse  amplitude  and  or  is  the  pulse  decay  constant.  Both 
of  these  terms  are  functions  of  altitude.  In  Ref.  4 it  is  shown  that 
the  noise  spectrum  produced  by  v such  pulses  occurring  each  second  is 
given  by 


P = A 


(?' 
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Further  experimental  work  described  In  Ref.  4 systematically 
explored  the  dependence  of  corona-dlscharge-nolse  characteristics  on 
altitude.  The  results  of  this  work  shown  In  Figure  3 demonstrate  that 
both  the  amplitude  and  spectral  character  of  the  noise  are  strongly 
altitude-dependent.  In  preparing  Figure  3,  curves  of  the  form  of  Eq.  (2) 
were  plotted  to  fit  the  experimental  data  as  well  as  possible. 

Previous  calculations  of  equivalent  noise  field  relied  on 
"look-up"  and  nomograph  techniques  to  obtain  P,  the  relative  source- 
spectrum  magnitude.  An  alternative  and  much  better  computer  solution  was 
found  by  defining  a,  v,  and  a 1,.  terms  of  an  analytical  fit  to  the  plotted 
dat:.,  which  was  a function  of  pressure. 

It  was  determined  that  a good  fit  to  A,  v,  and  a could  be  ob- 
tained by  expressing 


5 

0.25 

A = 7.90569  X 10 

P 

(3) 

3 

0.48 

V = 3.83767  X 10 

P 

(4) 

-2 

a = 2.7777  X 10 

P 

(5) 

where  p Is  In  torr. 

10 

A fit  to  the  1962  U.S.  Standard  Atmosphere  resulted  In  the 

expression 


p(torr) 


760  exp 


(h  -t-  0.002  h^\ 

, / 


(6) 


where  h is  in  kilofeet.  Inserting  Eq.  (6)  Into  Eqs.  (3),  (4),  and  (5) 
and  inserting  these  in  turn  into  tq.  (2)  re."=ulted  In  an  expression  for 
the  relative  source-spectrum  amplitude  that  was  uniquely  determined  by 
the  altitude  and  frequency  and  is  denoted  by  the  variable  name  PREL  in 
the  program. 
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RELATIVE  NOISE-CURRENT  SPECTRAL  DENSITY 


2, 


Streamer  Discharge 


In  the  program  of  Refs.  6 and  7,  flight-test  and  laboratory 
measurements  were  made  of  streamer  discharge  characteristics,  Tn  the 
work  of  Ref.  7,  an  analysis  was  made  of  the  signals  that  streamers  in- 
duce in  receiving  systems.  Certain  e.ssentials  of  this  analytical  work 
are  presented  in  Appendix  A along  with  extensions  and  generalizations 
required  for  the  present  program. 

The  results  of  the  streamer  studies  indicate  that  streamer 
discharges  induce  pulses  of  the  sort  shown  in  Figure  4,  which  illustrates 
the  pulses  generated  in  a wire  located  immediately  below  the  dielectric 
surface  on  which  streamers  occur.  It  is  evident  from  Figure  4 that  the 
time  wavefom  of  the  streamer  pulse  depends  on  the  length  of  the  streamer. 
The  amplitude  of  the  induced  pulse,  of  course,  depends  on  the  coupling 
between  the  streamer  and  the  receiving  "antenna." 

To  define  the  noise  spectrum  generated  by  streamers  occurring 
in  a practical  situation,  it  is  necessary  to  derive  a statistical 
description  for  a series  of  pulses  such  as  those  shown  in  Figure  4 oc- 
curring at  random  times  with  randcsnly  varying  amplitudes  and  time 
structures.  Such  statistical  analyses  were  carried  out  in  Appendix  A 
for  the  three  general  arrangements  of  dielect.’ic  region  and  receiving 
antenna  shown  in  Figure  A-1.  The  results  of  the  statistical  analyses 
are  given  in  Table  A-1  and  appear  to  be  quite  dissimilar.  However,  when 
the  statistical  results  are  used  to  calculate  noise-source  spectra  as 
shown  in  Figure  A-5,  it  is  found  that  the  resulting  noise-source  spectra 
generated  on  the  various  surfaces  are  virtually  Indistinguishable. 

In  a further  exploration  of  streamer-noise-source  character- 
istics, source-spectrum  calculations  were  carried  out  for  different 
sizes  of  the  dielectric  streamer-noise-source  region  (streamer  discharge 
lengths  varied).  The  results  of  these  calculations  indicate  that 
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FIGURE  4 TYPICAL  CURRENT  PULSES  INDUCED  BY  STREAMER  DISCHARGES 

Increasing  the  size  of  the  region  has  the  effect  of  increasing  the  low- 
frequency  content  of  the  streamer-noise-source  spectrum  as  shown  in 
Figures  A-6  and  A-7. 

Essentially  the  results  given  in  Appendix  A give  expressions 
for  streamer-noise-source  spectra  as  a function  of  charging  current  for 
a number  of  representative  geometrical  arrangements  of  dielectric  and 
antenna.  Of  particular  interest  in  the  results  of  Appendix  A are  the 
facts  that  source  spectral  characteristics  depend  primarily  on  the 
dimension  of  the  longest  streamer  that  can  occur  on  the  dielectric  region, 
and  that  they  do  not  depiend  critically  on  the  shape  of  the  dielectric 
region,  'ihis  means  that  in  the  computer  model  it  is  not  necessary  to 
be  precise  about  specifying  the  form  of  the  region  on  which  the  streamer 
discharges  occur. 
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C . Electrcmiagnetlc  Coupling 
1 . General 

Simply  specifying  the  spectral  character  of  the  noise  source 
as  was  done  in  the  previous  sections  is;  not  sufficient  to  determine  the 
noise  experienced  by  a system.  It  is  also  necessary  to  specify  the 
electromagnetic  coupling  between  the  noise  source  and  the  "antenna"  of 
the  affected  system.  Analytical  deterninations  of  coupling  in  simple 
physical  structures  were  made  in  Ref.  8.  An  experimental  technique  for 
measuring  coupling  was  developed  in  the  program  of  Ref.  7,  and  used  to 
study  coupling  between  a streamer  source  and  a nearby  antenna.  (In  this 
case  the  dimensions  of  the  coupling  structure  were  small  compared  to  the 
wavelength  of  signals  of  Interest.)  In  the  work  of  Ref.  4 the  method 
was  extended  to  the  determination  of  the  coupling  between  aircraft  ex- 
tremities where  corona  discharges  occur,  and  receiving-antenna  locations 
on  the  aircraft.  (Here  the  dimensions  of  the  structure  were  up  to  several 
wavelengths  long  at  the  frequencies  of  interest.)  Further  measurements 
were  made^^  to  determine  the  coupling  between  antennas  mounted  on  the 
fuselage  of  an  aircraft  and  streamers  occurring  on  the  radome  or  canopy. 

The  coupling -measurement  technique  has  involved  the  use  of  an 

isolated  broadband  spark  noise  souxce  positioned  at  the  p-statlc  source 

location,  and  either  an  oscillo.scope  or  a radio  receiver  connected  to 

the  antenna  in  question.  This  measurement  technique  yields  the  magnitude 

of  the  coupling  factor,  but  provides  no  information  regarding  phase. 

Using  the  coupling  data,  it  is  possible  to  specify  the  noise  spectrum  at 

an  antenna  locatlm,  but  g \erally  it  is  not  possible  to  predict  the 

time  waveform  of  the  noise  pulses  induced  in  the  antenna  terminals.  To 

calculate  time  waveforms,  it  is  necessary  to  know  both  the  phase  and 

12 

aiiiplitude  of  the  coupling  function,  (In  some  special  cases  it  is 


possible  to  use  supplemental  information  and  to  synthesize  an  expression 
for  the  phase  function.  In  general,  however,  the  available  coupling 
data  do  not  allow  pulse  waveform  prediction.) 

It  should  be  noted  at  this  time  that  the  systematic  developmei.t 
of  computer  techniques  such  as  the  method  of  moments  for  the  solution 
of  electromagnetic  problems  offers  the  promise  of  providing  information 
regarding  both  amplitude  and  phase  of  the  coupling  function.  To  in- 
vestigate the  utility  of  such  techniques,  a computer  run  was  made  using 
a program  already  available  at  SRI  to  calculate  coupling  fields  at  the 
extremities  of  a very  crude  cylinder  model  of  a swept-wing  aircraft. 

The  results  of  this  calculation  presented  in  Appendix  B are  very  en- 
couraging, since  the  computer  calculations  are  remarkably  accurate  in 
predicting  the  resonances  and  general  behavior  of  the  coupling  functions. 
Since  no  effort  was  made  to  match  the  size  of  the  "antenna"  used  in  the 
cOTiputer  calculations  with  that  of  an  actual  aircraft  antenna,  the  mag- 
nitude of  the  calculated  coupling  could  not  be  ccmipared  with  measured 
values.  There  is  no  reason  to  expect  that  the  magnitude  cannot  be 
predicted  with  accuracy. 

2 . KC-135  Corona  Source  Coupling 

a.  Measured  Coupling  Data 

The  physical  arrangement  considered  in  the  KC-135  noise- 
coupling measurements  of  Ref.  4 is  shown  in  Figure  5.  The  antennas  used 
in  making  the  measurements  were  a small  tail-cap  and  a flush  belly 
antenna  located  in  a fairing  at  the  root  of  the  wing.  Coupling  was 
measured  between  i?ach  of  these  antennas  and  the  noise  source  regions 
at  each  of  the  airfoil  extremities.  The  results  of  the  measurements 
are  shown  in  Figure  6. 
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FIGURE  5 ANTENNA  AND  CORONA-NOISE-SOURCE  LOCATIONS  FOR  KC-135  COUPLING 
MEASUREMENTS 


Although  the  coupling -measurement  technique  used  inevitably 
includes  the  test-antenna  characteristics  i,i  the  measured  results,  the 
form  of  the  coupling  as  a function  of  frequency  will  not  be  affected  by 
the  details  of  the  test  antenna,  provided  the  antenna  dimensions  are 
small  compared  to  a wavelength  at  the  frequencies  of  interest.  In  this 
case,  only  the  magnitude  of  the  coupling  function  will  be  affected  by 
a change  in  iiie  antenna.  The  data  shown  in  Figure  6 have  therefore  been 


adjusted  to  represent  the  coupling  to  an  antenna  having  unity  induction 
area,  a,  in  response  to  a low-frequency,  vertically  polarized  signal. 


The  variation  with  frequency  of  the  coupling  to  the 
various  points  is  of  considerable  interest,  since  it  shows  the  effect 
of  the  various  electromagnetic  lesonances  of  the  aircraft.  For  example, 
the  peak  in  coupling  between  the  belly  antenna  and  the  tips  of  the  tall 
surfaces  at  approximately  3 MHz  occurs  at  the  frequency  ■'hat  makes  the 
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FIGURE  6 MEASURED  COUPLING  FACTORS  FOR  ANTENNAS  ON  BOEING  367-80 
AIRCRAFT 


path  distance  from  a point  Just  aft  of  the  wings  to  a tip  of  one  of  the 

tail  surfaces  one-quarter  wavelength.  The  presence  of  the  resonance 

peaks  in  the  coupling  function  serves  to  emphasize  the  fact  that  the 

pulses  arriving  at  the  antenna  terminals  will  differ  from  the  pulses 

generated  by  the  corona  discharges  at  the  extremities  in  that  they  will 

12 

be  stretched  out  in  time  bv  aircraft  resonances. 


Also  of  Interest  is  the  variation  of  the  coupling  at  a 
particular  frequency,  with  position  along  an  airfoil.  This  is  illustrated 
in  Figure  7 which  is  a plot  of  the  variation  in  coupling  of  th>'  belly 
antenna  to  the  trailing  edge  of  the  wing  near  the  tip.  The  coupling  is 
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FIGURE  7 VARIATION,  NEAR  WING-TIP,  OF  COUPLING  TO 
BELLY  ANTENNA  IN  KC-135  AIRCRAFT 

proportional  to  the  charge  density  that  would  be  excited  by  the  antennas 
if  a voltage  were  applied  to  them.  At  all  frequencies  the  charge  tends 
to  concentrate  at  the  extremities,  as  shoA^n  by  the  plot  in  Figure  7. 
Thus,  the  coupling  tends  to  be  maximum  at  the  points  where  discharges 
are  most  likely  to  occur.  The  dip  in  coupling  at  11  MHz  occurs  at  a 
distance  approximately  one-quarter  wavelength  from  the  tip.  This  result 
also  Indicates  the  fallacy  in  assuming  that  the  noise  induced  in  an 
antenna  on  an  aircraft  will  always  decrease  as  the  noise  source  is 
moved  away  from  the  antenna.  From  Figure  7 it  is  apparent  that  on  the 
trailing  edge  of  the  wing  near  the  tip.  moving  a noise  source  away  from 
the  belly  antenna  inc.  ,ases  the  noise  induced  in  the  antenna, 

b.  Effects  of  Antenna  Size  on  Corona  Noise  Coupling 

The  data  of  Figure  6 can  be  applied  without  restriction 
to  aircraft  of  similar  shape  but  different  size.  In  general,  the  same 
resonances  will  occur  on  the  scaled  aircraft,  but  they  will  occur  at 
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the  scaled  frequenc’>^.  In  Ref.  4 it  is  shown  that  on  the  scaled  aircraft 

the  coupling  ^ (at  the  scaled  frequency  uid)  varies  as  the  square  root 
1/n 

of  the  scale  n of  the  aircraft  (n  > 1).  Thus, 


1/n 


t /n  . 


(7) 


(The  square  root  occurs  as  the  result  of  an  artifice  in  the  way  in  which 
the  absolute  corora-noise-source  level  is  defined.)  In  order  for  Eq.  (7) 
to  be  valid  everything  about  the  general  form  of  the  aircraft  (including 
antenna  size)  must  be  scaled.  This  means  that  the  induction  area,  a, 
of  the  test  antenna  varies  as  the  inverse  square  of  the  scale  n,  so  that 


a 


1/n 


(8) 


Thus  the  normalized  coupling  on  the  1/n-scale  aircraft  varies  as 

, 5/2 

Va  I = n Va  . (9) 

I 1/n 

3.  Coupling  to  Streamer  Sources 

The  spark  noise-source  technique  was  used  in  the  work  of  Ref.  11 
to  measure  the  coupling  between  three  antennas  located  on  the  fuselage 
of  an  aircraft,  and  the  surface  of  a radome  located  on  the  nose.  The 
measurements  were  made  using  an  oscilloscope  connected  to  the  antenna 
terminals,  and  represent  the  low-frequency  coupling,  which  is  of  primary 
interest  because  streamer  noise  spectra  fall  off  rapidly  with  increasing 
frequency  above  1 MHz.  The  results  of  these  measurements  are  reproduced 
in  Figure  8. 
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FIGURE  8 NOISE  COUPLING  TO  RADOME 

One  of  the  most  striking  features  evident  in  the  coupling  data 
of  Figure  8 is  the  fact  that  the  coupling  is  relatively  independent  of 
antenna  location  on  the  fuselage.  This  result  is  in  keeping  with  the 
work  in  Chapter  VII  of  Hef.  1 where  it  is  shown  that,  for  the  case  of 
a cylindrical  structure  with  a noise  source  at  one  end,  the  electric 
noise  field  intensity  under  quasi  static  conditions  (noito  coupling)  is 
independent  of  position  along  the  cylinder  see  Eq.  (61),  Ref.  <11.  Also, 
the  work  of  Ref.  1 indicates  that  the  magnitude  of  the  electric  field 
is  inversely  proportional  to  the  peripheral  distance  around  the  cylinder. 
This  relationship  is  useful  in  applying  the  data  of  Figure  8 to  aircraft 
of  different  size.  It  may  also  be  used  to  treat  antennas  located  on 
a tapered  section  of  the  aircraft  such  as  the  nose. 
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occurring.  First,  consideration  was  , ^ ven  to  the  idealized  geometry  of 
Figure  10,  which  also  allows  analytical,  -olution,  and  which  therefore 
also  permitted  evaluation  of  the  spark-discharge-probe  coupling- 
measurement  technique.  This  configuration  may  be  considered  representative 
of  the  practical  case  of  an  antenna  wire  mounted  on  the  back  of  a 
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dielectric  ..indov  cut  in  the  skin  ol  an  aircraft.  The  me''siir-cd  and 
calculated  values  of  couplinR  for  this  antenna  conf  ir ura  t ion  are  al.so 
shown  in  FiRure  10.  It  shvjuld  be  noted  that  the  coiiilinR  is  presented 
as  the  couplinR  factor  v.  not  normalized  by  ib<'  ...ii^'iia  induction  area 
as  was  done  earlier.  PresentinR  couplir.R  in  terms  of  ■ s appropriate 
in  this  case  since  the  current  induced  in  the  12-RaRi'  wiri'  by  a sti'eamer 
discharRe  will  not  be  aflected  by  the  leuRth  ol  tlie  wire.  This  is  true 
because  the  rcRion  of  the  antenna  wire  in  the  immediate  vicinity  of  tlu' 
streamer  i , most  important  ii;  coupliiiR  streamer  siRiials.  It  should 
not  be  inferred  from  this  that  the  streamer  noise  is  indi'pcndent  ol  tlu' 
IcnRth  of  the  antenna  and  dielectric  rcRion.  A lotiRer  rcRion  means  a 
higher  total  charging  current  and  consequently  a higher  prl  (.)  and  ;■ 
higher  streamer  noise  level. 
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FIGURE  10  CALCULATED  VALUES  OF  COUPLING  FACTOR  COMPARED  TO  VALUES 
MEASURED  BY  FIELD  PROBE 

Although  a compietf?ly  idealized  geometrical  arrangement  such 
as  that  illustrated  in  Figure  10  is  not  likely  to  be  found  on  a flight 
vehicle,  the  structure  is  sufficiently  sL-nilar  to  the  case  of  a slot  cut 
in  a vehicle's  skin  that  the  results  are  of  Interest,  It  !'■  evident 
that  the  coupling  is  quite  uniform,  with  a value  of  roughly  3 m ^ over 
most  of  the  gap,  and  th  t it  increases  only  in  the  high-coupling-f ield 
regions  in  the  vicinity  of  the  wire  and  the  edge  of  the  sheet. 

Also  considered  in  Kef,  7 was  the  coupling  geometry  of  one 
of  the  special  radome  antennas  used  in  the  B-47  flight  tasks  in  that 
program  and  illustrated  in  Figure  A-l(a).  Here,  because  of  the  three- 
dimensional  nature  of  the  geometry,  it  was  found  that  the  coupli'g  varies 
witti  position  on  the  radome  surface.  It  was  found,  however,  that  the 
magnitude  of  the  coupling  is  of  the  same  order  as  in  the  case  of  the 
idealized  geometry  of  Figure  10.  Indeed,  the  average  streamer  coupling 
for  the  radome  cf  Figure-  ,\-l(a)  is  roughlv  f - 3m 
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D.  Equivalent  Noise  Field 


Although  the  terminal  noise  current  completely  specifies  the  pre- 
cipitation-static noise  response  of  an  antenna,  terminal  current  is  not 
a famili£r  term.  Furthermore,  since  the  terminal  current  induced  in  an 
antenna  depends  on  the  antenna  characteristics,  the  terminal  current  is 
of  no  general  value  unless  the  antenna  characteristics  are  also  specified. 
Thus  the  precipitation-static  noise  data  would  be  more  convenient  if  they 
were  expressed  in  terms  of  the  field  intensity  required  to  generate  the 
calculated  terminal  current  in  the  antenna.  Although  the  antenna  noise 
current  depends  on  size  and  sensitivity,  the  equivalent  noise  field 
applies  to  any  antenna  installed  at  the  specified  location.  In  particular, 
one  generally  knows  the  desired  signal  field-strength  environment  in 
which  his  system  must  function.  If  he  also  is  able  to  express  p-static 
noise  in  terms  of  equivalent-noise-f leld  strength,  he  can  immediately 
specify  the  ratio  of  signal  to  precipitation-static  noise  in  his  -stem. 

Consider  first  the  low-frequency  case  in  which  a dipole-type  antenna 

mounted  on  the  aircraft  will  have  the  radiation  pattern  of  an  elementary 
13 

dipole.  The  characteristics  of  the  antenna  are  detemiined  by  specify- 
ing the  direction  of  the  dipole  axis  and  the  induction  area,  a,  of  the 
antenna,  which  is  defined  as 


a 


(10) 


where  E is 
1 

the  induced 
Bolljahn  in 
antenna  are 
capaci tance 


the  field  intensity  of  the  incident  wave,  q is 

-9 

charge  on  the  antenna,  and  = (1/36:)  x 10 
Ref.  13  shows  that  h (effective  height)  and  a 


related  through  the  expression  e a - h C whei 

0 e 

of  the  antenna. 


the  value  of 
farad/meter, 
for  a given 
C is  the 
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For  the  low-frequency  case  being  considered,  the  current  flowing 

at  the  short-circuited  antenna  is  jx<!.  Thus,  using  Eq.  (10)  we  may  write 

I = j'JUG  • Rearranging  and  making  an  analogy  between  noise  and 
sc  0 1 

signal,  we  obtain 


(11) 


where  I is  the  tenuinal  noise  current,  and  E is  defined  as  the  equlv- 
n n 

alent  noise  field. 

As  the  frequency  is  raised  above  the  quasi-static  range,  the  antenna 

can  no  longer  be  represented  by  a simple  dipole.  In  calculating  |e  I, 

I nl 

therefore,  it  is  necessary  to  account  for  changes  in  the  radiation 
pattern  and  for  the  effects  of  aircraft  resonances.  In  Ref.  4,  equivalent 
noise  fields  at  the  KC-135  fin  cap  for  frequencies  of  2 MHz  and  higher 
were  calculated  using  probe-antenna  noise  currents  together  with  probe- 
antenna  data.  Since  antenna  gain  varies  with  direction  of  signal  arrival, 
the  equivalent  roise  field  calculated  in  Ref.  4 did  not  have  a unique 
value  for  a given  frequency,  but  was  a function  of  direction.  A com- 
parison of  the  results  of  noise-field  calculations  using  radiation 
patterns  with  results  from  quasi-static  theory  indicated  that  the  error 
stemming  from  the  use  of  the  simpler  quasi-static  theory  is  not  greater 
than  that  resulting  from  small  changes  in  direction  of  propagation. 

Thus,  it  was  concluded  in  Ref.  4 that  the  quasi-static  method  provides 
a simple  method  for  estimating  the  order  of  magnitude  of  equivalent 
noise  field  at  a given  antenna  location  when  pattern  and  impedance  data 
for  an  antenna  at  that  location  are  not  available.  Quasi-static  es- 
timates of  antenna  properties  are  used  in  this  reprart  in  calculating 
noise  fields. 
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E.  Aircraft  Charging 


1.  Total  Aircraft  Charging  Current  (Corona) 

To  properly  proceed  with  the  planning  of  a discharger  Instal- 
lation, It  Is  necessary  to  estimate  the  charging  that  one  should  expect 
on  the  aircraft  under  Investigation.  Since  extensive  flight-test  and 
laboratory  studies  of  charging  were  made  on  the  707  prototype.  It  Is 
convenient  to  structure  the  analysis  so  that  these  existing  data  can  be 
Interpolated  and  applied  to  the  aircraft  In  question.  Actually,  the 
approach  Is  to  estimate  the  ratio  of  the  new  aircraft  current  to  the 
707  current  when  the  two  aircraft  ere  operated  In  the  same  precipitation 
conditions.  To  provide  a concrete  Illustration,  the  procedure  will  be 
carried  out  here  for  the  case  of  the  F-4  aircraft. 

The  precipitation  charging  current  to  a vehicle  Is  given  by 

1 = q c V A (12) 

p eff 


where 

q = Charge  per  particle 
P 

c = Particle  concentration  In  cloud 
V = Aircraft  velocity 

A = Effective  Intercepting  area  of  aircraft, 
eff 

Since  we  have  assumed  that  the  two  aircraft  are  operating  In 
the  same  precipitation,  the  first  three  terms  on  the  right  of  Eq.  (12) 
will  be  the  same  for  the  two  aircraft  when  they  are  operated  at  the 
same  speed.  In  this  case,  the  ratio  of  the  charging  currents  will  be 
given  simply  by  the  ratio  of  the  effective  Intercepting  areas. 
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If  the  aircraft  operate  at  different  speeds,  the  situation  Is 

n 

only  slightly  more  complicated.  Laboratory  tests  confirmed  by  flight 
14 

tests  Indicate  that  the  Individual  particle  charge  q varies  only 

P 

slightly  with  speed  In  the  speed  range  of  present  day  aircraft.  A 

elf 

will  also  vary  somewhat  with  speed. 

Let  us  now  consider  the  question  of  determining  the  effective 

Intercepting  area  of  an  aircraft  of  Interest.  Although  considerable 

analytical  work  has  been  done  to  predict  Intercepting  areas  of  aircraft 

flying  In  clouds  composed  of  spherical  water  droplets,  the  ability  to 

obtain  these  solutions  Is  critically  dependent  on  the  fact  that  the 

drag  coefficient  of  a moving  sphere  can  be  determined.  These  results 

cannot  be  applied  directly  to  the  case  of  an  aircraft  operating  In  a 

cloud  COTiposed  of  Irregularly  shaped  Ice  crystals  for  which  the  drag 

coefficient  cannot  be  determined.  Thus,  In  estimating  Ice-crystal 

Impingement  on  aircraft.  It  Is  necessary  to  rely  heavily  on  flight-test 

experimental  results  and  to  use  the  results  of  water-droplet  impingement 

calculations  to  Indicate  trends  as  various  parameters  such  as  size  are 

varied.  The  procedure  is  not  satisfying.  In  that  considerable  fitting 

and  approximating  must  be  done.  (Until  more  analysis  or  flight  testing 

is  carried  out,  no  other  approach  Is  possible.)  For  example,  it  Is 

indicated  on  page  54  of  Ref.  4 that  the  effective  area  of  the  707  Is  a 

small  fraction  of  the  projected  area  (roughly  20%  to  25%  Interception 

efficiency).  On  page  119  of  the  same  report.  It  Is  observed  that  the 

percent  of  the  frontal  area  of  an  ellipsoid  effective  In  intercepting 

15-17 

particles  varies  inversely  with  the  length  of  the  body.  From  this 

it  was  concluded  in  Ref.  4 that  in  the  case  of  similar  aircraft,  the 
ratio  of  the  effective  areas  A is  given  by  the  ratio  of  the  spans  of 
the  aircraft  S: 
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(13) 


A S 

F-4  _ F-4 

A ■ S 
707  707 

Later  the  problem  of  scaling  was  reviewed,  and  it  was  noted 
that  the  frontal  area  of  the  airfoils  was  generally  almost  equal  to  the 
frontal  area  of  the  fuselage.  It  was  argued,  then,  that  since  the  per- 
cent effectiveness  of  the  airfoils  is  much  higher  than  the  effectiveness 
of  the  fuselage,  the  charging  would  stem  primarily  from  the  airfoils. 

Thus,  the  way  in  which  aircraft  charging  is  affected  by  changes  in  size 
is  determined  by  the  way  in  which  airfoil  interception  area  varies  with 

size.  Considerably  after  the  publication  of  Ref.  4 the  authors  discovered 
18 

NACA  RM  No.  9A05,  in  which  a study  of  the  impingement  of  water  droplets 
on  several  airfoils  was  reported.  The  data  in  this  NACA  report  (repro- 
duced in  Figure  16)  indicate  that,  in  the  vicinity  of  20%  to  25%  inter- 
ception efficiency,  percent  intercepting  area  varies  inversely  with 
airfoil  scale.  As  indicated  earlier,  flight-test  data  obtained  on  the 

Boeing  367-80  aircraft  (the  prototype  of  the  707  and  KC-135)  indicated 

2 

that  during  much  of  the  time  the  intercepting  area  was  50  ft  (25%.  of 
the  airfoil  frontal  area)  or  less  as  shown  in  Table  III  of  Ref.  4. 

These  last  observations  confirm  that  the  relationship  of  Eq.  (13)  is 
correct  for  estimating  the  effects  of  scalir"  on  aircraft  charging  current. 
Careful  Inspection  of  Figure  11  indicates  that  the  percent  intercepting 
area  varies  inversely  with  aircraft  scale  up  to  roughly  a 257o  inter- 
ception efficiency,  and  that,  for  small  aircraft,  the  percent  inter- 
ception grows  iiore  slowly  with  decreasing  size.  Thus  for  the  F-4 
operating  at  KC-135  speeds,  the  charging  current  would  be  somewhat  lower 
than  predicted  by  Eq,  (13).  However,  since  the  F-4  can  operate  at 
supersonic  speeds,  the  appropriate  free-stream  Reynold's  number  for  the 
F-4  case  in  Figure  11  will  be  higher  than  for  the  KC-135  (twice  as  high 
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for  M = 2) , This  has  the  effect  of  continuing  the  straight-line  portion 
of  the  curves  in  Figure  11  to  a scale  modulus  appropriate  for  the  F-4. 
Accordingly,  it  is  appropriate  to  use  Eq.  (13)  in  estimating  F-4  charg- 
ing current  from  XC-135  flight-test  data. 

The  wing-spans  of  the  F-4  and  KC-135  are  as  follows: 

S = 3S.3  ft 

F-4 

S = 13t  ft 

707 

Using  these  values  in  Eq.  (13)  we  obtain 

■f-4  = '707  • 

Values  of  total  charging  current  calculated  in 
in  the  estimation  of  corona-noise  levels  to  be 
and  in  devising  discharge  installations. 

2 . Localized  Charging  (Streamers) 

To  determine  the  noise  levels  generated  by  streamer  discharges, 
it  is  necessary  to  be  able  to  specify  the  charging  current  responsible 
for  generating  the  streaming.  Attempting  to  make  precise  calculations 
of  dielectric  charging  current  is  even  more  difficult  than  calculating 
total  charging  current.  Accordingly,  we  will  adopt  the  procedure  of 
using  empirical  estimates  based  on  flight-test  measurements  made  on  the 
F-4  tests  described  in  Ref.  14.  The  test  aircraft  Included  provisions 
for  measuring  the  current  arriving  on  the  windshield  and  on  isolated 
patches  on  the  radome. 

In  formulating  the  computer  model,  for  the  windshield  and 
canopy  we  will  use: 


(14) 

this  manner  are  required 
expected  on  an  airframe 
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, front,  diel. 

I = i I 

diel.  A T aircraft 

front,  aircraft 


(15) 


where  I is  the  current  arriving  on  the  windshields  or  canopy, 

diel. 

^ ,,  , is  the  frontal  area  of  the  dielectric  region  in  question, 

front,  diel. 

A„  , „ is  the  projected  frontal  area  of  the  aircraft,  and 

front,  aircraft 

I . „ is  the  total  charging  current  to  the  aircraft. 

T aircraft 


For  the  radrane  we  will  use: 


forward  3 ft 


radome  A T aircraft 

front,  aircraft 


(16) 


where  I is  the  current  arriving  on  the 

radome 

is  the  projected  frontal  area  of  the  forward 


radome  and  A 

forward  3 ft 
three  feet  of  the  radome. 
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Ill  COMPUTTT?  MO!^ELING  OF  EXISTING  DATA 


A . General 

The  experimental  and  analytical  data  regarding  p-static  noise  dis- 
cussed in  Section  II  were  processed  appropriately,  and  made  part  of  the 
computer  program.  In  some  cases,  this  meant  using  an  analytical  expres- 
sion or  approximation  to  the  desired  parameter.  In  other  cases  where 
the  data  did  not  lend  itself  to  approximation,  it  was  simply  stored  in 
tabular  form. 

A brief  description  will  be  given  here  of  the  approach  taken  in 
order  that  the  reader  be  familiar  with  the  contents  of  the  program  and 
how  it  ./as  derived. 

B.  Corona  Noise 

The  total  charging  current  arriving  on  the  aircraft  is  calculated 

using  Eqs.  (12)  and  (13),  where  q and  c have  been  replaced  by  an  expres- 

P 

Sion  relating  th /se  quantities  to  aircraft  speed  and  cloud  type  (e.g., 

cirrus,  s tratacumul us , frontal  snow).  The  free  paiameters  in  this 

expression  have  been  adjusted  such  that  the  charging  current  to  a 707 

flying  in  the  same  frontal  snow  conditions,  I , is  4 mA.  This  adjust- 

707 

ment  was  made  to  correspond  to  the  highest  charging  conditions  encounter- 
ed on  the  707  prototype  during  the  flight- test  program  of  Ref.  4.  The 
statistical  data  presented  in  Figure  27  of  Ref.  4 indicate  that  this 
charging  is  exceeded  only  C.0027o  of  the  flight  time  on  a jet  transport. 
Thus,  the  noise  levels  obtained  from  the  computer  program  (whi -h  are 
related  co  the  charging  current)  will  be  typical  of  the  most  severe 
chaxging  conditions  for  a given  altitude,  speed,  and  type  of  precipitation 

Precedinf  paga  blank 


particle.  The  program  also  includes  a printout  of  the  probability 
that  the  actual  charging  current  will  exceed  the  calculated  charging 
current . 

Provisions  are  included  in  the  program  to  allocate  the  discharge 
currents  among  the  various  aircraft  extremities.  Based  on  the  work  of 
Ref.  4 the  currents  are  apportioned  approximately  as  follows; 

Each  empennage  tip  1/7 
Each  wing-tip  2/7 

Spccifyinjj  "tiis  discliargG  i*Gglons  in  coinputot*  codG  was  n©cGSsary  sine© 
the  coupling  to  the  noise  source  depends  critically  on  the  aircraft 
extremity  from  which  the  corona  discharge  is  occurring. 

Provision  is  also  made  for  quieting  selected  discharge  locations. 
This  quieting,  essentially  a 40-dB  noise  reduction,  allows  the  electro- 
magnetic compatibility  (EMC)  engineer  the  opportunity  to  study  the 
effect  of  installing  corona  quenching  or  quieting  devices  (p-static 
discharges  for  example)  at  selected  locations. 

Once  the  corona  current  discharged  from  an  extremity  is  determined, 
Eqs.  (2)  through  (6)  are  used  to  calculate  the  source  spectrum  of  the 
RF  noise  that  this  discharge  generates.  Tlie  received  noise  spectrum 
resulting  from  the  discharge  at  the  extremity  in  question  is  calculated 
by  convolving  the  noise-source  spectrum  with  the  appropriate  coupling 
function  from  Figure  6 for  the  antenna  location  and  extremity  under 
consideration.  The  total  received  noise  spectrum  is  obtained  by  summing 
the  contributions  from  each  of  the  extremities  on  an  rms  basis. 

Since  the  coupling  data  of  Figure  6 are  experimental  and  do  not 
lend  themselves  to  an  easy  analytical  descripi ion , they  are  read  into 
the  program  from  data  cards  and  stored  in  the  program  in  tabular  form. 
This  technique,  while  somewhat  cumbersome,  has  the  advantage  that  wlien 
additional  coupling  data  become  available,  they  can  be  used  and  no  pro- 
gram modifications  are  necessary. 
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The  program  includes  provisions  to  scale  the  frequencies  of  the 
resonances  in  the  coupling  curves  of  Figure  6 to  make  them  appropriate 
for  the  aircraft  under  investigation.  Eq.  (9)  is  used  to  adjust  the 
magnitude  of  the  coupling  function  according  to  the  size  of  the  aircraft. 

Once  the  totr.l  received  noise  spectrum  is  calculated,  Eq.  (11)  is 

used  to  express  the  received  noise  signal  in  terms  of  an  equivalent 

noise  field  at  the  antenna  location.  It  should  be  observed  at  this 

time  that  corona-noise-coupling  data  are  available  only  for  an  aircraft 

the  shape  of  a KC-135.  However,  the  same  data  have  been  applied  to 

14 

aircraft  such  as  the  F-4,  which  is  quite  dissimilar  from  the  KC-135, 
with  remarkably  good  results.  Thus,  until  additional  coupling  data 
are  generated,  there  is  no  choice  but  to  use  existing  data.  As 
described  above,  provisions  arc  included  in  the  program  to  incorporate 
new  coupling  information  as  it  oecomes  available. 


C.  Streamer  Noise 

The  current  arriving  on  the  dielectric  surface  under  consideration 

is  calculated  using  Eq.  (15)  for  the  ease  of  a structure  such  as  a 

windshield  or  canopy,  and  using  E(|.  (16)  loi’  a radoine.  To  use  these 

equations,  it  is  necessary  to  know  the  total  aircraft  cliarging  current 

I . This  is  available  from  the  earlier  corona-noise  calculaticns. 

T 

Noise-source  calculal i('ns  are  made  using  Eq,  (A-3)  of  Appendix  A. 

In  the  first  term  of  this  equation,  I and  v are  simply  constants 

m 

related  to  the  characteristics  of  a streamer  discharge,  and  the  PKF  ' 
is  found  from  the  dielectric  charging  current  by 

■ = ■ <’7, 

1.5  X 10 
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The  second  term  of  Eq.  (A-3)  is  given  by  Eq,  (A-5) , in  which  all  the 
terms  except  ui  are  constants  related  to  the  characteristics  of  the 
streamer  discharge.  The  last  term  in  Eq,  (A-3)  is  given  in  Table  A-1 
for  three  different  geometric  arrangements.  Since  it  was  shown  that 
the  noise-spectrum  characteristics  are  not  sensitive  to  the  shape  of 
the  dielectric  region,  the  computer  program  uses  the  third  expression 
in  Table  A-1.  In  this  case  the  term  W is  the  length  in  me;ers  of  the 
minimum  characteristic  dimension  of  the  dielectric  region  under  consid- 
eration— i.e.,  the  width  of  a rectangular  region,  (Thi  term  W may  be 
thought  of  as  roughly  twice  the  length  of  the  louge >t  possible  streamer 
discharge  on  the  region  in  question.) 

It  should  be  noted  that  the  last  expression  in  Table  A-1  includes 
the  coupling  factor  which  must  be  specified.  Three  different  physi- 
cal arrangements  of  "antenna"  and  streamer  source  can  be  specified. 

The  first  is  that  illustrated  in  Figure  8,  where  the  antenna  is  located 
on  the  fuselage  aft  of  the  region  (such  as  the  radome)  In  which  the 
streamering  is  occurring.  The  data  of  Figure  8 have  been  stored  in  the 
program,  so  that  it  is  necessary  only  to  specify  the  diameter  F of  the 
fuselage  at  the  antenna  location  and  the  distance  d the  noise  source  is 
located  ahead  of  the  metallic  portion  of  the  fuselage. 

The  second  streamer  configuration  is  illustrated  in  Figure  9,  where 
the  antenna  is  located  on  the  fuselage  beneatli  a windshield  or  canopy 
on  which  streamers  are  occurring.  It  is  evident  from  Figure  9 that  the 
coupling  is  relatively  independent  of  the  location  d of  the  discharge 
on  the  canopy.  Accordingly,  tho  data  stored  in  the  program  use  the 
approximation  that  the  coupling  is  entirely  independent  of  d.  It  is 
necessary  only  to  specify  the  distance  C of  llie  antenna  aft  of  tlie  front 
edge  of  the  canopy, 

!t  should  be  noted  that  in  both  Figure  8 and  Figure  9,  the  noise 
source  is  sufficiently  far  removed  from  the  antenna  that  the  noise 
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coupling  will  be  affected  by  antenna  induction  area  a.  Thus  the  noise 

coupling  is  presented  in  normalized  form  ('i'/a)  and  it  is  possible  to  use 

Eq.  (11)  directly  to  calculate  at  equivalent  noise  field  E . 

n 

The  third  streamer  coupling  configuration  considered  is  illustrated 
in  Figure  10.  Here  a wire  is  located  on  the  underside  of  the  dielectric 
on  which  streamers  are  occurring.  With  this  arrangement  the  streamer 
noise  source  is  so  close  to  the  wire  "antenna"  that  the  noise  picked  up 
by  the  region  of  the  wire  in  the  immediate  vicinity  of  the  streamer 
dominates.  Thus  the  streamer  coupling  is  independent  of  wire  length 
(i.e.,  independent  of  "antenna"  induction  area  a).  In  this  case  the 
computer  out.vjt  indicates  the  short-circuit  current  induced  in  the  wire 
lather  than  the  equivalent  noise  field  at  the  wire  location. 

D.  Verification  of  Computer  Model 

1 . General 

To  test  the  functioning  and  accuracy  of  the  computer  modeling, 
calculations  were  carried  out  to  predict  precipitation-static  noise 
levels  under  various  conditions.  Vehicle  types  and  sizes  were  chosen 
to  correspond  as  closely  as  possible  to  the  vehicles  employed  in  the 
flight  tests  described  in  Refs.  4 through  7,  and  14.  Precipitation 
charging  conditions  also  were  chosen  to  duplicate  as  nearly  as  possible 
the  conditions  encountered  during  the  flight  tests.  During  the  flight 
tests,  measurements  were  made  of  both  streamer  and  corona  noise,  and 
three  different  types  of  aircraft  were  used.  Thus,  it  is  possible  to 
verify  the  computer  model  under  a gratifyingly  large  variety  of  con- 
ditions. 
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2. 


Corona  Noise 


a . 707  Belly  Antenna 

In  the  flight  tests  of  Refs,  4 and  5,  the  707  prototype 
aircraft  was  operated  without  any  p-static  noise-reduction  provisions  in 
light  cii  us  c.'ouds  such  that  the  total  charging  current  to  the  aircraft 
was  250  uA.  Measurements  were  made  of  noise  current  induced  in  the 
belly  antenna.  These  measured  data  have  been  interpreted  as  equivalent 
noise  field,  and  are  plotted  as  the  five  data  points  in  Figure  12, 


FIGURE  12  CORONA-NOISE  FIELDS  AT  707  BELLY  ANTENNA 


Also  plotted  in  Figure  12  is  the  curve  of  equivalent 
noise  field  at  the  707  belly  antenna  location  predicted  bv  the  computer 
model.  In  the  computing  of  the  predietotl  noise  field,  the  program  was 
artifically  constrained  to  a charging  curient  of  250  -A  typical  of  the 
very  light  cirrus  clouds  encountered  during  this  portion  of  the 
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flight- test  program.  (The  computer  model  would  normally  estimate  a 
charging  current  of  1 mA  to  a 707  aircraft  in  cirrus.) 

It  is  evident  from  Figure  12  that  the  computer  model  is 
remarkably  accurate  in  predicting  all  the  details  of  the  corona  noise 
induced  in  the  belly  antenna.  The  absolute  level  of  the  noise-field 
intensity  is  accurately  predicted,  and  the  airframe  resonances  are 
correctly  represented. 

b.  707  Tail-Cap  Antenna 

In  the  707  flight-test  program  of  Refs.  4 and  5,  noise 
measurements  were  also  made  on  the  tail-cap  antenna  when  the  aircraft 
was  operated  in  light  cirrus  such  that  the  total  aircraft  charging 
current  was  250  uA.  In  the  first  flight  tests,  the  aircraft  was  entirely 
unmodified  so  that  corona  discharges  could  occur  from  each  of  the  air- 
foil extremities.  The  noise-field  intensities  measured  in  these  flights 
are  shown  as  circles  in  Figure  13.  Also  shown  in  the  upper  part  of  the 
figure  are  the  noise  fields  at  the  tail  cap  predicted  by  the  computer 
model.  (As  before,  the  computer  program  was  constrained  to  a charging 
current  of  250  uA  typical  of  the  extremely  light  cirrus  clouds  encountered 
during  this  test  series.)  It  is  evident  from  the  upper  part  of  Figure  13 
that  the  computer  code  accurately  predicted  the  magnitude  and  spectral 
character  of  the  noise  in  the  tail-cap  antenna  when  all  corona  discharges 
were  permitted. 

Later  in  the  same  series  of  flight  tests,  electrical 
insulating  tape  was  applied  to  the  trailing  edge  of  the  rudder  of  the 
707  test  aircraft.  In  this  way  corona  discharges  from  the  rudder  were 
inhibited.  The  noise  fields  measured  at  t.ie  tail-cap  antenna  under 
these  conditions  are  shown  as  the  dots  in  Figure  13.  Shown  for  com- 
parison in  the  lower  pai’t  of  the  figure  is  a curve  giving  the  noise 
spectrum  predicted  by  the  computer  model  for  tail-cap  antenna  noise 
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FIGURE  13  CORONA-NOISE  FIELDS  AT  707  TAIL-CAP  ANTENNA 


when  discharges  from  the  rudder  are  inhibited.  Again  it  is  evident 
that  the  computer  code  accurately  predicted  both  the  magnitude  and  the 
character  of  the  noise  spectrum. 

The  utility  of  being  able  to  request  that  the  computer 
model  suppress  one  of  the  noise  sources  is  evident  in  Figure  13.  Here 
it  is  evident  that  noise  generated  by  discharges  from  the  rudder  domi- 
nates the  tail-cap  antenna  noise  environment.  By  eliminating  the 
rudder  corona,  the  noise  level  is  reduced  by  almost  40  dB  at  the  lower 
frequencies.  Such  information  is  highly  significant  to  the  EMC  engineer 
since  it  means  that  in  developing  noise-reduction  techniques,  great 
care  must  be  exercised  to  make  certain  that  rudder  corona  noise  is 
eliminated . 

In  comparing  the  data  of  Figure  12  with  the  lower  curve 
in  Figure  13  we  see  that  the  noise  levels  are  comparable,  and  that  the 
affec  '.s  of  airfram  ? resonances  are  evident  in  the  peaks  of  noise  levels 
at  frequencies  of  roughly  3 MHz  in  both  sets  of  data.  It  should  also 
be  observed  that,  in  these  two  cases,  the  noise  levels  are  substantial, 
and  that  functioning  of  radio  systems  could  be  improved  by  further  noise 
reductions. 

To  further  demonstrate  the  utility  of  the  computer  model, 
a calculation  of  tail-cap  antenna  noise  was  carried  out  for  the  case  of 
the  707  aircraft  operating  in  frontal  snow  where  the  aircraft  charging 
current  is  4 mA,  The  resulting  noise  levels  shown  in  Figure  14  are 
extremely  severe — even  with  rudder  discharges  inhibited. 

c.  F-4  Tail-Cap  Antenna 

A further  test  of  the  computer  model  was  its  ability  to 
predict  noise  levels  on  an  aircraft  much  smaller  than  the  707.  Accord- 
ingly, the  computer  code  was  used  to  predict  the  noise  levels  to  be 
expected  at  the  tail  cap  of  an  F-4  aircraft  operating  in  cirrus  clouds, 
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EQUIVALENT  NOrSE  FIELD  IN  1-lcHz  BANDWIDTH — /iV/m 


FIGURE  14  PREDICTED  CORONA-DISCHARGE  NOISE  FIELDS  ON  707  TAIL-CAP 
ANTENNA  IN  FRONTAL  SNOW 


It  was  recognized  that  the  shape  of  an  F-4  is  greatly  different  from 
that  of  a 707,  and  that  coupling  data  gathered  using  a model  of  the 
707  was  not  appropriate  for  the  F-4,  so  that  good  agreement  between 
measured  and  predicted  data  should  not  be  expected.  The  computer 
results  are  shown  in  Figure  15  along  with  data  measured  in  the  F-4 
flight-test  program  of  Ref.  14. 

It  is  noteworthy  in  Figure  15  that  eliminating  corona 
discharges  from  the  rudder  affords  considerably  less  noise  reduction 
than  it  did  on  the  707  (see  Figure  13).  This  stems  from  the  fact  that 
the  F-4  is  sufficiently  small  that  coupling  fran  the  fin-cap  to  the 
wing  and  elevator  tips  is  almost  as  high  as  the  coupling  from  the  tail- 
cap  to  the  rudder.  Thus  rudder  corona  noise  is  not  as  dominant  a 
source  of  interference  on  the  F-4. 

It  is  also  evident  that  the  degree  of  agreement  between 
predicted  and  measured  noise  in  Figure  15  is  not  as  good  as  it  was  for 
the  707.  This  is  reasonable  since  the  coupling  data  used  in  the  present 
computer  model  are  not  strictly  appropriate  for  the  F-4.  Instead  of 
concentrating  on  the  differences  between  the  predicted  and  measured 
results  in  Figure  15  it  is  more  important  to  observe  that  the  predicted 
data  are  still  in  remarkably  good  agreement  with  the  measurements.  In 
fact  the  difference  between  predicted  and  measured  results  is  of  the 
same  order  as  the  difference  caused  by  flying  the  aircraft  in  a 
different  cloud  type — e.g.  , stratociimulus  ’.  stead  of  cirrus. 

In  general,  we  should  conclude  that  the  computer  model 
does  a remarkably  good  job  of  predicting  the  noise  levels  on  an  aircraft 
differing  considerably  from  the  707  in  both  shape  and  size. 
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FIGURE  15  CORONA-NOISE  FIELDS  AT  F-4  TAIL-CAP  ANTENNA 


d.  F-4  Belly  Antenna 

As  an  additional  test  of  the  general  applicability  of 
the  computer  model,  it  was  used  to  predict  the  noise  level  on  the  belly 
of  the  F-4  aircraft.  The  resul  s of  this  calculation  are  shown  in 
Figure  16.  It  can  be  argued  that,  because  of  the  general  symmetry  of 
the  aircraft's  structure,  the  belly  data  should  also  be  applicable  to 
the  case  of  an  antenna  mounted  on  the  top  of  the  fuselage.  Accordingly, 
data  measured  using  a Chelton  "towel-bar"  antenna  mounted  on  the  top 
of  the  fuselage  of  the  F-4  of  Ref.  14  are  shown  in  Figure  16  for  com- 
parison. Again,  it  is  evident  that  the  computer  output  provides  an 
acceptable  estimate  of  the  noise  actually  measured  in  flight. 

Thus,  we  should  conclude  from  Figure  16  that  the  computer 
model  can  be  stretched  considerably  and  applied  to  situations  otner  than 
those  precise  ones  for  which  there  are  measured  coupling  data. 

3.  Streamer  Noise 

a.  B-47  Coated  Radome  Antenna 

The  streamer-noise  portion  of  the  computer  model  was 
checked  by  comparing  its  predictions  with  flight-test  data  generated  on 
the  B-47  aircraft  used  in  tlie  tests  of  Refs.  6 and  7.  The  B-47  was 
equipped  with  a s ecial  nose  section  divided  into  separate  metal  cavities 
with  an  acrylic  radotn.'  section  forming  the  outside  wall  of  each  of  the 
cavities.  Two  of  these  cavities,  located  on  opposite  sides  of  the 
fuselage,  were  eacu  fitted  with  a short  antenna  wire  which  measured  the 
noise  at  its  location.  The  acrylic  cover  over  the  first  of  these  cavity 
sections  was  coated  with  electrically  conducting  paint  to  prevent 
streamering  on  its  surface.  The  second  cavity  cover  was  left  uncoated 
so  that  electric  charge  could  accumulate  on  the  acrylic  surface,  thereby 
causing  streamer  discharges  to  occur  on  the  uncoated  radome  section. 
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EQUIVALENT  NOISE  FIELD  IN  1-kHz  BANDWIDTH  — /JV/m 


FIGURE  16  CORONA-NOISE  FIELDS  AT  F-4  CHELTON  ANTENNA  (top  of  fuselage) 
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The  streamer  noise  generated  on  the  uncoated  section  coupled  into  its 
antenna  and  also  into  the  antenna  under  the  coated  section  on  the  oppo- 
site side  of  the  aircraft. 

The  coupled  streamer  noise  current  measured  at  the  coated 
radome  section  is  shown  in  Figure  17. 
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FIGURE  17  NOISE  SPECTRUM  IN  COATED  RADOME  SECTION  ON  B-47  AIRCRAFT 
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Also  shown  in  the  figure  are  predicted  noise  currents  generated  by  the 
p-static  computer  model  for  the  case  of  an  antenna  mounted  on  an  air- 
craft such  that  F = 1 ft  and  d = 2 inches  for  the  case  of  cirrus  and 
stratocumulus  clouds.  It  is  evident  that  the  predicted  noise  currents 
are  in  good  agreement  with  the  measured  data. 

b.  B-47  Uncoated  Radome  Antenna 

Noise  measurements  were  also  made  on  the  antenna  under 
the  uncoated  radome  section.  This  arrangement  corresponds  to  the  case 
of  an  antenna  such  as  those  often  mounted  on  the  inside  surface  of  a 
pilot's  canopy.  Here  the  antenna  wire  is  positioned  on  the  back  surface 
of  the  frontal  dielectric  on  which  streamers  are  occurring.  Obviously 
the  coupling  in  this  condition  is  high,  so  that  the  induced  noise 
currents  should  be  high. 

Noise  currents  measured  in  the  antenna  behind  the  uncoated 
radome  surface  are  shown  in  Figure  18,  Unfortunately,  as  is  indicated 
in  Ref.  7,  following  the  flight  tests  it  was  found  that  the  preamplifier 
used  with  this  antenna  could  be  overdriven  by  the  individual  streamer 
noise  pulses  so  that  the  measured  noise  spectrum  is  not  truly  representa- 
tive of  that  existing  at  the  antenna.  In  particular  the  low-frequency 
content  of  the  measured  noise  was  demonstrated  to  be  accentuated.  In 
spite  of  this  shortcoming  in  the  measured  data,  it  is  interesting  to 
note  that  the  measured  noise  currents  in  the  uncoated  antenna  are 
several  orders  of  magnitude  higher  than  those  measured  under  the  coated 
radome . 

Also  plotted  in  Figure  18  are  the  short-circuit  noise 
currents  predicted  by  the  computer  model  for  the  uncoated  radome 
antenna.  The  general  noise  levels  are  in  good  agreement  with  the 
measured  data.  The  disparity  in  the  shapes  of  the  measured  and 
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FIGURE  18  NOISE  SPECTRUM  IN  UNCOATED  RADOME  ANTENNA  SECTION  ON  B-47 
AIRCRAFT 


computed  spectra  should  undoubtedly  be  attributed  to  the  shortcomings 
of  the  B-47  flight-test  instrumentation. 

It  should  be  concluded,  therefore,  that  the  computer 
model  is  able  to  predict  streamer  noise  spectra  with  good  accuracy,  and 
that  the  predicted  results  can  be  applied  with  confidence. 


IV  APPLICATION  AND  LIMITATIONS 
OF  COMPUTER  MODEL 


A.  General 

Since  the  computer  model  developed  here  is  entirely  based  on  the 
results  of  earlier  work,  ito  applications  and  limitations  are  circom- 
scribed  by  this  earlier  work.  In  each  of  the  previous  programs,  an 
effort  was  made  to  ensure  that  the  results  were  as  generally  applicable 
as  possible.  However,  practical  considerations  dictated  that  the  flight 
tests  and  supporting  laboratory  work  restrict  themselves  to  the  con- 
sideration of  a limited  number  of  aircraft.  As  a result  the  existing 
noise-source  data  are  of  general  applicability,  since  they  are  indepen- 
dent of  vehicle  design.  Corona  coupling,  on  the  other  hand,  has  been 
measured  only  for  the  KC-135  aircraft  (which  was  used  as  the  flight- 
test  vehicle).  The  KC-135  corona  coupling  data  can  be  scaled  with  no 
loss  of  accuracy  to  similar  aircraft  of  different  size,  and  can  be 
applied  with  reduced  accuracy  to  aircraft  of  other  shapes.  For  corona 
noise,  supporting  flight-test  data  are  available  for  the  KC-135  and  for 
the  F-4  aircraft. 

The  streamer-noise  data,  while  also  somewhat  limitad,  are  of  quite 
general  applicability  since  they  do  not  depend  on  aircraft  design  and 
RF  resonance  to  the  extent  that  the  corona  noise  coupling  does. 

B.  Appl ications 

1.  Corona  Noise 

Although  corona-noise  coupling  measurements  were  made  only  for 
the  case  of  a fin-cap  located  near  the  rear  of  the  vertical  tail  and 
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for  an  AUF  sense  antenna  or  the  belly  of  the  KC-135,  the  results  can  be 
applied  somewhat  more  generally.  For  example,  the  fin-cap  coupling  data 
can  be  applied  with  good  accuracy  to  an  antenna  located  anywhere  on  the 
top  of  the  vertical  fin.  These  might  range  from  a probe  antenna  such 
as  is  used  on  the  707  to  an  electrically  isolated  portion  of  the  fin 
structure  such  as  that  used  on  the  UC-8. 

The  belly-antenna  coupling  data  were  measured  for  the  case  of 
an  antenna  located  at  the  wing  root  as  Indicated  in  Figure  5.  In  view 
of  the  symmetry  of  the  aircraft,  however,  these  data  also  can  be  applied 
with  little  loss  of  accuracy  to  the  case  of  an  antenna  located  at  the 
same  fuselage  station  on  the  top  of  the  fuselage.  If  necessary,  the 
same  coupling  data  can  be  used  with  reasonable  accuracy  to  investigate 
coupling  to  antennas  located  as  much  as  one-quarter  fuselage  length 
fore  or  aft  of  the  wing-root  position. 

Although  no  coupling  measurements  were  made  for  the  case  of  an 
antenna  located  on  a wing-tip,  it  is  possible  to  do  a reasonable  job  of 
predicting  noise  levels  on  these  too.  Looking  at  Figure  6 we  note  that 
for  the  tail-cap  antenna,  the  highest  noise  coupling  is  to  discharges 
occurring  from  the  rudder.  We  may  argue  with  reasonable  validity  that 
the  case  of  a wing- tip  antenna  the  dominant  source  of  corona  noise 
will  be  corona  discharges  from  the  wing-tip  trailing  edge.  Since  a 
wing-tip  is  indistinguishable  electrically  from  a fin-tip,  the  coupling 
from  wing-tip  antenna  to  wing-tips  should  be  given  quite  accurately  bv 
the  rudder-tip- to- tail-cap  coupling  curve  of  Figure  6. 

A similar  sort  of  argument  is  possible  for  the  case  of  noise 
coupled  into  a wing-tip  antenna  from  discharges  on  the  rudder.  Here 
we  can  invoke  reciprocity  and  argue  that  the  coupling  is  given  by  the 
wing-tip-to-tai 1-cap  coupling  of  Figure  6. 
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In  the  case  of  noise  generated  on  the  opposite  wing-tip  and  on 
the  elevator  tips,  unfortunately,  we  cannot  make  use  of  the  existing 
coupling  data  to  determine  coupling  factors. 

Since  only  the  coupling  data  of  Figure  6 exist  at  present,  they 
must  be  applied  to  all  aircraft  regardless  of  shape  or  size.  In  the 
case  of  extremes  such  as  a delta-wing  fighter,  fur  exan.ple,  it  would 
seem  that  the  validity  of  using  the  KC-135  coupling  data  is  doubtful. 

It  must  be  recalled,  however,  that  these  data  were  reasonably  accurate 
in  predicting  the  corona  noise  levels  on  an  F-4  aircraft  in  Ref.  14. 

In  using  the  program  to  study  corona  noise,  it  should  be 
recalled  that  one  of  the  outputs  is  the  equivalent  noise  field  at  the 
"antenna"  location.  This  information  can  be  used  directly  as  an  input 
to  the  Introsystem  EMC  Analysis  Program  to  analyze  the  noise  coupled 
into  the  "antenna"  (a  wire,  an  aperture,  or  even  a real  intentional 
a.Ttenna)  at  the  location  in  question. 

In  the  case  of  situations  involving  radio  receiving  systems, 
the  equivalent  noise  field  at  the  antenna  location  is  often  the  only 
parameter  of  interest.  One  can  immediately  determine  whether  the 
resulting  signal-tc'-noife  ratio  is  acceptable.  If  the  p-static  noise 
is  unacceptably  high,  a knowledge  of  the  noise-field  intensity  indicates 
to  the  EMC  engineer  the  degree  of  corona  noise  reduction  required  by 
his  system. 

2.  Streamer  Noise 

On  most  aircraft  of  standard  design,  the  fraction  of  tlie 
frontal  surfac,-’  fabricated  of  highly  insulating  dielectric  material  has 
been  small.  This  does  not  mean  that  the  frontal  surfaces  have  always 
been  fabricated  of  metal.  It  simply  means  that  where  dielectric  materi- 
als were  used,  it  was  almost  always  possible  to  coat  the  dielectric  with 
a conductive  layer  such  as  a metal  film  or  a conductive  paint.  Thus 
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the  total  current  involved  in  streamer  discharges  has  been  much  smaller 
than  the  corona  current,  and  consequently  the  streamer-noise-source 
intensity  is  generally  proportionally  lower  than  the  corona-noise-source 
intensity.  Accordingly,  it  has  been  possible  to  restrict  consideration 
of  streamer  noise  to  sources  in  the  close  vicinity  of  the  receiving 
"antenna".  For  this  reason  available  streamer  coupling  data  have  been 
generated  exclusively  for  the  case  of  antennas  located  near  the  di- 
electric siiriace  on  which  the  streamer  dischages  are  occurring.  For- 
tunately, constraints  on  aircraft  design  are  such  that  these  data  are 
adequate  to  cover  most  cases  of  practical  interest. 

The  main  frontal  regions  of  dielectric  on  an  aircraft  are 
the  radome,  the  windshield,  and  the  canopy  structure.  The  PSTAT  pro- 
gram is  designed  to  treat  the  coupling  of  noise  from  streamers  on  these 
surfaces  to  "antennas"  located  on  the  fuselage  of  the  aircraft.  The 
output  of  the  program  is  designed  to  provide  either  the  equivalent 
noise  field  at  the  location  of  interest  or  the  short-circuit  noise 
current  induced  in  the  structure  of  interest,  (To  obtain  the  short- 
circuit  current  it  is  necessary  to  specify  the  "antenna"  induction  area. 
These  outputs  permit  the  investigation  of  streamer  effects  on  intention- 
al antennas  of  receiving  systems  as  well  as  noise  coupling  into  inadver- 
tent antennas  formea  by  wires,  cables,  etc. 

In  some  cases  the  "antenna"  wire  is  located  immediately  below 
the  surface  on  which  streamers  are  occurring.  Provisions  are  included 
in  the  PSTAT  program  to  treat  these  areas.  With  this  geometry,  the 
signal  induced  in  the  wire  is  determined  primarily  by  coupling  to  the 
region  in  the  immediate  vicinity  of  the  streamer.  Therefore  the  noise 
pulses  are  largely  independent  of  antenna  length  so  that  a calculation 
of  equivalent  noise  field  is  not  appropriate.  Consequently,  the  out- 
put  of  the  program  provides  the  short-circuit  noise  current  induced  in 
a wire  on  the  back  side  of  a dielectric.  The  wire,  of  course,  may  be 
an  antenna  wire,  or  cabling  associated  with  an  aircraft  avionic  system. 
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C.  Limitations 


1 . Corona  Noise 

Since  the  computer  program  v/as  based  on  an  existing  technolo- 
gical base,  it  can  be  no  better  than  the  avilable  base.  Accordingly, 
there  are  restrictions  on  the  regions  of  the  aircraft  for  which  noise- 
level  calculations  can  be  made.  For  the  same  reason,  the  coupling  data 
are  most  accurate  for  aircraft  shaped  like  the  KC-135.  No  coupling 
data  exist  for  vehicles  such  as  rockets. 

Since  the  coupling  data  were  measured  with  a system  that  did 
not  provide  phase  information,  only  the  magnitude  of  the  coupling  factor 
is  included  in  the  program.  For  this  reason,  the  program  provides  only 
noise-power  spectral-density  information.  It  is  not  possible  at  present 
to  calculate  amplitudes  of  individual  noise  pulses.  (Such  noise-pulse 
information  is  important  in  treating  digital  systems.)  To  calculate 
individual  noise  pulses  it  is  necessary  to  have  information  on  both  the 
phase  and  magnitude  of  the  coupling  function. 

Coupling  information  were  generated  only  to  frequencies  of 
12  MHz  in  the  KC-135  measurement  program.  Thin  sets  the  upper  limit  on 
the  frequency  range  that  may  be  considered. 

All  of  the  corona-noise  spectral  data  depend  on  pulse  shape 
and  spectral  measurements  made  on  the  program  of  Ref.  4.  These  measure- 
ments were  carried  out  with  state-of-the-art  instrumentation  available 
at  that  time.  More  recently,  preliminary  investigations  using  more 
modern  spectrum-analyzer  techniques  appear  to  indicate  that  the  spectrum 
does  not  decay  as  rapidly  in  the  VHF  region  as  was  indicated  by  the 
work  of  Ref.  !.  Over  the  region  where  corona  noise  is  most  severe,  the 
noise-source  data  appear  to  be  accurate. 
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2.  Streamer  Noise 

As  in  the  case  of  corona  noise,  the  computer  model  was  founded 
on  the  existing  technological  base.  The  geometries  for  which  coupling 
data  exist  allow  consideration  of  a reasonably  wide  variety  of  problems. 

It  is  likely,  however,  that  situations  may  arise  for  which  coupling 
data  do  not  exist.  Fortunately,  the  problem  here  is  not  as  complex  as 
in  the  case  of  corona-noise  coupling,  since  the  st’’eamer  source  must 
generally  be  close  to  the  affected  "antenna"  so  that  vehicle  resonances 
are  not  important  in  determining  the  coupling. 

The  streamer  noise-source  data  are  based  on  oscillographic 
measurements  of  waveforms  made  during  the  programs  of  Refs.  6 and  7, 

It  is  very  likely  that  there  is  more  high-frequency  structure  associated 
with  the  propagation  of  a streamer  than  was  evident  on  the  instrumenta- 
tion available  at  that  time.  The  spectrum  out  to  frequencies  of  roughly 
10  MHz  is  adequately  represented  by  the  existing  data. 
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V CONCLUSIONS  AND  RECOMMENDATIONS 


The  objectives  of  the  present  program  were  to  employ  existing 
information  regarding  precipitation-static  noise  generation  and  coupling 
processes  to  develop  a computer  model  that  could  be  incorporated  into 
the  Intrasystem  EMC  analysis  program.  Such  a model  was  generated,  and 
tested  by  comparing  its  predictions  with  flight-test  data  generated  on 
earlier  p-static  research  programs.  Within  the  limits  of  existing 
data,  it  was  possible  to  predict  with  good  accuracy  the  precipitation- 
static  noise  levels  to  be  expected  in  flight. 

The  computer  program  includes  provisions  for  suppressing  the  noise 
from  one  or  more  of  the  various  sources  to  permit  the  EMC  engineer  to 
investigate  the  importance  of  applying  noise-reduction  measures  to 
each  of  the  sources.  This  capability  is  important  because  implementing 
noise-reduction  techniques  on  certain  parts  of  the  aircraft  may  be 
costly,  whereas  in  other  areas  they  can  be  implemented  with  ease.  Thus, 
the  EMC  engineer  is  always  in  the  position  of  having  to  defend  the 
necessity  lor  carrying  out  all  of  his  noise-reduction  recommendations. 

In  preparing  the  computer  model,  every  effort  was  made  to  incorpo- 
rate the  maximum  flexibility  commensurate  with  avoiding  undue  complexity. 
An  effort  was  made  to  recognize  areas  where  growth  of  the  model  should 
be  expected,  and  to  provide  for  the  possibility  of  future  expansion  in 
these  areas. 

To  use  the  resulting  computer  model,  the  EMC  engineer  must  prepare 
a set  of  input  cards  describing  the  aircraft,  and  its  speed,  altitude 
regime,  location  r>f  receiving  "antennas",  etc.  The  output  of  the  pi’ogram 
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provides  information  regarding  the  total  charging  current  to  the  air- 
craft together  with  the  magnitudes  of  the  noise  fields  at  the  locations 
of  interest  as  well  as  the  magnitudes  of  the  short-circuit  currents 
induced  in  antennas  of  specified  induction  area  at  these  locations.  In 
general,  these  outputs  are  adequate  to  specify  the  noise  environment 
sufficiently  to  permit  judgements  regarding  the  functioning  of  avionic 
systems  and  regarding  the  need  for  fixes. 

As  was  recognized  at  the  outset,  the  computer  model  is  limited  by 
the  available  technological  base.  This  became  more  apparent  as  develop- 
ment of  the  computer  model  progressed.  In  particular  there  is  a need 
for  more  information  on  noise  coupling.  It  is  important  to  be  able  to 
specify  the  coupling  to  more  than  two  general  regions  on  a particular 
aircraft.  Of  even  greater  importance,  the  coupling  data  need  to  be 
extended  to  other  vehicle  types.  These  should  include  a wide  variety 
of  aircraft  types  as  well  as  rockets,  helicopters,  etc. 

The  increasing  Importance  of  digital  avionics  systems  indicates 
that  attention  should  be  given  to  their  susceptibility  to  individual 
noise  pulses.  The  present  coupling  data  do  not  contain  information 
about  the  phase  of  the  coupling  function,  so  that  it  is  not  possible 
to  predict  the  structure  of  individual  p-static  noise  pu?.ses  coupled 
into  systems.  In  future  work  to  generate  coupling  data,  the  scheme 
used  should  include  the  capability  of  providing  phase  information. 

The  coupling  data  base  can  be  enlarged  either  by  experimental 
measurements  or  by  computational  means.  Recent  developments  in  both 
instrumentation  a id  computer  methods  offer  considerable  promise  in  this 
regard.  Ideally,  one  should  provide  for  both  analytical  as  well  as 
experimental  work,  to  generate  faith  in  the  results. 

On  the  present  program  an  effort  was  made  to  investigate  the 
feasibility  of  using  computational  techniques  such  as  the  method  of 
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moments  to  generate  coupling  data.  The  very  crude  model  used  was 
remarkably  accurate  in  predicting  the  form  and  magnitude  of  the  coupling 
between  the  tail-cap  antenna  region  and  a noise  source  on  an  elevator 
tip  on  the  KC-135  aircraft.  These  techniques  should  definitely  be 
exploited  D'urther,  since  they  provide  phase  as  well  as  amplitude  infor- 
mation as  a matter  of  course. 

The  process  of  preparing  the  computer  model  served  to  reemphasize 
the  fact  that  the  basic  p-static  noise-source  data  are  based  on  measure- 
ments made  up  to  20  years  ago,  and  are  consequently  limited  by  the 
instrumentation  available  at  that  time  as  well  as  by  the  investi^oator' s 
view  of  what  was  required.  In  particular,  limits  on  available  instrumen- 
tation made  it  difficult  to  define  the  character  of  the  noise  sources 
above  10  to  20  MHz.  This  was  not  considered  a limitation  at  that  time 
because  VHF  systems  were  operated  at  relatively  high  signal  levels. 

Also,  microwave  system  design  was  such  that  it  was  quite  impervious  to 
p-static  interference.  At  present,  with  satellite  communication  systems 
being  considered  for  aircraft,  it  will  be  necessary  for  VHP’  systems  to 
operate  at  substantially  lower  signal  levels,  so  that  VHF  ccmiponents  of 
p-static  noise  may  be  significant.  For  these  reasons,  the  noise-source 
data  should  be  examined  using  modern  instrumentation. 

In  view  of  the  increased  use  of  dielectric  materials  in  aircrnft 
construction,  a careful  review  should  be  made  of  the  streamer  noise 
problem.  The  spectral  data  in  this  area  are  based  on  the  most  primitive 
of  tie  p-static  flight-test  measurements. 

A further  consequence  of  the  increased  use  of  composite  materials 
is  that  the  electromagnetic  shielding  provided  by  metal  skin  will  no 
longer  exist.  This  result  is  immediately  recognized  in  the  case  cf 
frontal  dielectric  surfaces  exposed  to  precipitation  charging  and  the 
consequent  streamering.  If  large  regions  of  dielectric  are  incorporated 
into  the  design,  however,  the  reduced  shielding  may  also  manifest  itself 


in  Increased  susceptibility  to  corone  noise.  The  consequences  of  such 
reduced  shielding  must  be  assessed. 

Whenever  possible,  the  future  work  should  include  provisions  for 
flight-test  verification  of  the  results.  Such  testing  is  essential  to 
provide  confidence  in  the  analytical  and  laboratory  experimental  results, 
as  well  as  an  opportunity  to  uncover  unexpected  effects.  Ideally,  the 
aircraft  selected  for  the  tests  should  include  as  many  as  possible  of 
the  materials  and  systems  recognized  to  lead  to  p-static  susceptibility. 
These  should  Include  composite  surfaces,  and,  particularly  avionics 
systems  employing  digital  subsystems  and  computers  operating  with  low- 
level  signals.  The  tests  should  include  measurements  of  noise  source 
characteristics  using  modem  instrumentation,  and  monitoring  of  the 
functioning  and  susceptibility  of  the  aircraft’s  avionics  systems. 
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Appendix  A 

STREAMER  DISCHARGE  CHARACTERISTICS 
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Appendix  A 


STREAMER  DISCHARGE  CHARACTERISTICS 


1.  General 

In  Kef.  7 the  results  of  laboratory  work,  flight  test  measurements 
on  a B-4  aircraft,  and  theoretical  analysis  are  combined  on  page  47  to 
yield  an  expression  for  the  short-circuit  current  induced  in  a receiving 
antenna  by  the  occurrence  of  a streamer  discharge  on  a plastic  surface. 
The  result  [Eq.  (26)  of  Ref.  7]  is 


where  u(z)  is  the  unit  step  defined  by 

u(z)  =0  for  z < 0 
u(z)  = 1 for  z >0, 

'i  defines  the  electromagnetic  cou -ling  between  the  streamer  source  and 
the  receiving  system,  and 

I , a,  b,  i,  B,  1,  ■' 
m 
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are  all  characteristics  of  the  streamer  breakdown  itself.  The  values 

I 

of  i and  V are  as  follows: 

I = Streamer  length 

5 

V = Streamer  velocity  = 1.27  X 10  m/s. 

and  the  remainder  are  obtained  from  the  expression  for  the  current  flow- 
ing past  a point  in  the  streamer  discharge  packet  [Eq,  (21)  of  Ref.  7], 

jT£pl*C<LiuC€d  il0Z*£  • 


{ 


II/  \ / -at  -l^t\ 

'2  (^’  *)  ■ Ux(‘  ^ * j 

where,  for  a typical  streamer, 

a = 0.597 
b = 0.403 

7 -1 

>=  1.67  X 10  s 

a 6 -1 

3 ^ 0.47  X 10  s 
I = 0.01  A. 


Equation  (A-2)  was  evaluated  for  a typical  value  of  coupling  factor, 

1'  3 m (typical  for  a wire  located  immediately  below  the  streamer), 

and  for  se.^eral  streamer  lengths,  and  the  results  plotted  in  Figure  4 
(Section  II-B).  Although  the  time  waveform  of  individual  noise  pulses 
induced  by  streamer  discharges  is  of  interest,  it  is  necessary  to  be  able 
to  express  the  results  in  terms  of  the  noise  spectrum  generated  by  a 
series  of  these  pulses  occurring  at  random  times  with  varying  discharge 
lengths  and  appropriately  differing  amplitudes.  Such  a calculation  was 
carried  out  in  Ref,  7,  and  numerical  results  were  obtained  for  the  case 
of  streamers  occurring  on  a triangular-shaped  radome  section  directly  ovt,' 
an  antenna  wire.  (In  this  case,  if  and  <*  both  varied  from  discharge  to 
discharge . ) 
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The  general  expression  for  the  short-circuit  current  noise-power 
spectrum  induced  in  an  antenna  by  a series  of  streamer  pulses  is  given 
by: 

2 2 

vl  u 2 

G(uu)  = |g(uu)|  |f(x,i)|2  . (A-3) 

where,  as  before,  I is  the  peak  value  of  the  current  flowing  in  the 
m 

discharge,  u is  the  velocity  of  propagation  of  the  discharge,  and  v 
is  the  average  rate  of  occurrence  of  the  streamers.  The  term  |g(u))|^  is 
the  square  of  the  Fourier  transform  of  the  basic  function 


a 

a 


(‘ 


(A-4) 


used  in  describing  the  pulse  produced  by  an  average  streamer  ^see 
Eq.  (A-1)]  and  is  given  by 


g(i) 


2 2 

1 I (a  + b)  + (a3  + bo-) 


(A-5) 


I |2 

where  the  parameters  are  defined  as  in  Eq.  (A-2).  The  term  |f(x,0| 
in  eq.  (A-3)  is  the  ensemble  average  of  the  function 


(A -6) 


over  the  random  variables 
discharges  occurring  on  a 


x and  i’  , and  defines  the  statistics  of  the 
particular  surface. 


To  perform  the  indicated  averaging,  it  will  be  necessary  to  know 
the  probability  density  h^(:',x)  of  a discharge  of  length  t originating 
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DIELECTRIC  REGION 
ON  WHICH 
CHARGING  OCCURS 

ANTENNA  WIRE 
ON  BACK  OF 
DIELECTRIC 


(a)  TRIANGULAR  DIELECTRIC,  ANTENNA  WIRE  ON  BACK 


DIELECTRIC  REGION 
ON  WHICH 
CHARGII4G  OCCURS 


RECEIVING 

ANTENNA 

DISTANCE 

>D 


(b)  TRIANGULAR  DIELECTRIC,  ANTENNA  DISTANT 


DIELECTRIC  REGION 
ON  WHICH 
CHARGING  OCCURS 


RECEIVING 

ANTENNA 

DISTANCE 

>W 


(c)  RECTANGULAR  DIELECTRIC,  ANTENNA  DISTANT 


SA-2919-22 


FIGURE  A-1  DIELECTRIC  SHAPES  AND  ANTENNA  ARRANGEMENTS  FOR  WHICH 
STREAMER  STATISTICAL  CALCULATIONS  WERE  CARRIED  OUT 
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frcm  the  point  x on  the  rim  around  the  charged  dielectric  region.  These 
statistics  may  be  inferred  frtxn  the  manner  in  which  the  dielectric  sur- 
face becomes  charged.  Such  an  evaluation  was  carried  out  in  Ref.  7 for 
a special  radcme  used  on  the  B-47  flight-test  aircraft.  The  calculations 
of  Ref.  7 will  be  reproduced  here  and  will  be  extended  to  other  gecmetries 
of  interest  so  that  ultimately  we  will  have  expressions  for  the  three 
combinations  of  dielectric  surface  and  antenna  location  illustrated  in 
Figure  A-1.  The  arrangements  shown  in  the  figure  were  chosen  because 
they  represent  the  range  of  charging  situations  and  antenna  arrangements 
one  is  likely  to  encounter  in  practice.  Figures  A-l(a)  and  (b)  are 
typical  of  the  charging  and  breakdown  situation  one  would  have  either 
on  a fighter  windshield  or  on  a radome  equipped  with  lightning-diverter 
strips.  The  shape  of  Figure  A-l(c)  corresponds  to  a dielectric  region 
located  on  the  leading  surface  of  an  airfoil. 

The  arrangement  of  Figure  A-l(a)  is  representative  of  an  "antenna" 
located  immediately  below  the  dielectric  surface  on  which  charging  and 
streamers  occur.  Figures  A-l(b)  and  (c)  arc  representative  of  situations 
in  which  the  antenna  element  is  not  immediately  adjacent  to  the  dielectric 
surface  on  which  streamering  is  occurring. 

2.  Evaluation  of  Streamer  Statistics 

a.  Triangular  Dielectric  Over  Antenna  Wire 

In  the  flight  tests  of  Ref.  7,  a special  nose  consisting  of 
four  triangular -shaped  radome  sections  was  installed  on  the  nose  of  the 
B— 17  flight-test  aircraft.  For  purposes  of  this  analysis  it  will  be 
convenient  to  approximate  the  radome  section  by  a plane  triangle  and  to 
assume  that  all  of  the  charging  of  the  radome  section  occurred  in  the 
unshaded  region  shown  in  Figure  A-2.  Another  assumption  that  will  be 
made  at  this  time  is  that  the  entire  unshaded  region  of  Figure  A-2 
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FIGUnE  A-2  APPROXIMATION  TO  REGION 
OF  RADOME  IN  WHICH 
CHARGING  OCCURRED 


charged  at  a uniform  rate.  Although  the  charging  rate  of  the  radome 
surface  did  var,  ith  distance,  x,  from  the  nose  of  the  radome,  the 
evaluation  of  th ' i variation  is  a lengthly  aerodynamic  problem.  Since 
very  little  is  kno\/n  regarding  the  aerodynamic  properties  of  ice  crystals, 
the  results  of  the  aerodyncamic  calculation  would  only  be  of  a qualitative 
nature.  Thus  it  was  felt  that  the  increased  accuracy  that  might  result 
from  including  the  variation  of  charging  rate  in  the  present  calculations 
was  not  worth  the  additional  effort. 

Since  it  has  been  assumed  that  the  unshaded  region  of  Figure  A-2 
charges  at  a uniform  rate,  the  probability  of  a discharge  occurring  from 
a point  0 < X < x written  P(0  < x)J  is  proportional  to  the  area 

a = 2 x(x  tan  (A-7) 

enclosed  by  the  triangle  MNP  shown  in  the  figure.  The  total  area  in 
whicii  charging  occurs  is  given  by 
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A = tan  0 


(AS) 


Thus , 


2 

P(0<X<x)=-  = — 0<x<D 

A 2 
D 

0 elsewhere. 

The  probability  density  h^  (x)  of  a discharge  from  a point  x along  the 
rim  is,  therefore, 


h3(x) 


^ _ 2x 

dx  2 
D 


= 0 


0 < X < D 


elsewhere. 


(A-9) 


Given  that  the  discharge  occurs  frcm  the  point  a discharge 
extending  to  any  point  on  the  line  PN  in  Figure  A-2  is  equally  likely 
since  it  has  been  assumed  that  the  entire  surface  charges  at  a uniform 
rate.  Thus  the  conditional  probability  density  of  a discharge  of  length 
Lwritten  h^(i’|x)]  is  given  by  the  rectangular  density  function 

, 1 

h ( x)  = 0 < f < X tan  6 

2 ' X tan  0 

(A-10) 

= 0 elsewhere. 

The  joint  density  function  h^(C,  x)  is  given  by 
h^(  r,  x)  -=  h^(£|x)h^(x) 
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which,  upon  substituting  Eqs.  (A-9)  and  (A-10) , becomes 


h^(i,  x)  = ~ 


D tan  0 


Letting  y = tan  9,  we  obtain 


h^(i,  X)  = — 
D V 


= 0 


0 < I < yx 
0 < X < D 

elsewhere. 


(A-11) 


With  the  above  Information  available,  we  may  now  evaluate  |f(x, 
Reproducing  Eq,  (A-6) , 


f(x,  1) 


- ♦(«) 


we  may  write 


f(x,  l)\^  = ,!,^(x) 


which  becomes 


(A-6) 


t(x,  ol^  = 2*^(x)  h - ros  --  jl 


/,  « \ 

V ” / 

The  ensemble  average  over  the  random  variables  x and  £ is  given  by 


(A-12) 


|f(x,  Z) 


/'  vx  D 

^-0  J x=0 


M (x,  O I h^(£,  x)dxdf 


Substituting  Eqs,  (A-12)  and  (A-11)  we  obtain 
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4 


|f(x,  l)\^ 


2 

D V 


X) 


dxd£  . 


The  integration  with  respect  to  i may  be  performed  immediately  to  give 


f(x,  l)\^ 


dx 


(A-13) 


It  will  now  be  necessary  to  obtain  an  expression  for  'li(x). 

The  measured  values  of  4 one  centimeter  in  from  the  rim  in  Ref.  7 are 

shown  plotted  versus  x in  Figure  A-3,  It  is  evident  frcm  the  figure 

that  the  points  lie  very  nearly  along  the  straight  line,  'ir  = mx  where 
-2 

m = 32  m . 


0 0.05  0.10  0.15  0.20 

DISTANCE  FROM  NOSE,  X — m 
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FIGURE  A-3  VARIATION  OF  COUPLING  FACTOR  WITH 
DISTANCE  FROM  RADOME  NOSE 


gives 


Substituting  the  above  expression  for  Kx)  into  Eq.  (A-13) 


(A-14) 
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Upon  performing  the  indicated  integration  and  rearranging  terms 
Eq.  (A-14)  becomes 


f(x,  l)\  = 


2 2 
2lD  » 

*■  <T" 

Y UJ 


VD 
sin — w 
2v 


yD 

2v 


0) 


yD 

Sin  0) 
YD  V 

+ COS  uu  - 2 

u yr* 

0) 

u 


(A-15) 


b.  Triangular  Dielectric,  Distant  Antenna 

Since  we  are  considering  a dielectric  surface  that  has  the 
same  shape  and  exposure  to  precipitation  as  in  Section  2-a,  above,  the 
arguments  regarding  the  statistics  of  streamer  discharge  location  and 
length  will  be  unchanged.  Thus,  we  may  simply  start  with  Eq.  (A-13) , and 
modify  the  succeeding  steps  as  appropriate  for  the  different  antenna 
coupling. 

Because  the  receiving  antenna  is  not  in  the  immediate  vicinity 
of  the  dielectric  surface  and  surrounding  rim,  the  coupling  will  be  far 
less  sensitive  to  position  than  it  was  in  the  previous  analysis.  Let  us 
assume,  in  fact,  that  the  coupling  is  constant  and  independent  of  x in 
Figure  A-2.  Substituting  into  Eq.  (A-13)  we  obtain 


f(x,  f)i 


dx.  (A-16) 


Upon  performing  the  Integration  and  rearranging,  Eq.  (A-16)  becomes 


f(x,  O = 2'if 


1 - - 


\D 

sin  ^ uu 
2u 


2v 


(A-17) 
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c. 


Rectangular  Dielectric,  Distant  Antenna 


Let  us  now  consider  the  case  of  a rectangular  dielectric  region 
as  illustrated  in  Figure  A-4  undergoing  i liform  charging.  The  probability 


1 

i 

IF) 

1 

1 

1 

0 
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FIGURE  A-4  RECTANGULAR  DIELECTRIC  REGION  OF  CHARGING 


of  a discharge  occurring  from  a point  0 < X < x [written  P(0  < X < x)J 
is  proportional  to  the  area 


a = Wx  (A-18) 

enclosed  by  the  rectangle  to  the  left  of  the  dotted  line  in  the  figure. 
The  total  area  in  which  charging  occurs  is  given  by 


Thus , 


A = LW  . 


(A-19) 


P(0  < X < x)  = - 
A 

= - 0 < X < L (A-20) 


= 0 elsewhere  . 


The  probability  density  h^(x)  of  a discharge  from  the  point  x along  the 
rim  is,  therefore, 
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h3(x) 


dP 

dx 


1_ 

L 


= 0 


0 < X < L 
elsewhere 


(A-21) 


Given  that  the  discharge  occurs  from  the  point  x,  a discharge  extending 
to  any  point  along  the  line  QR  is  equally  likely,  since  it  has  been 
assumed  that  the  entire  surface  charges  at  a uniform  rate.  Thus,  the 
conditional  probability  density  of  a discharge  of  length  I [written 
h^C-Cjx)]  is  given  by  the  rectangular  density  function 


h^C  «|x) 


W 


= 0 


0<£<- 
elsewhere  , 


(A-22) 


The  joint  density  function  h^(£,  x)  is  given  by 

h^(£,  x)  = h^ii’lx)  h^Cx) 

_ i 

~ W L (A-23) 

2 W 

= — for  0 < ^ < — 

LW  2 

0 X < L 


Picking  up  the  argument  at  Eq.  (A-12) 

|f(x,  l)\^=  2v^(x)  (1  - cos  7 1)  . (A-12) 

The  ensemble  average  over  the  random  variables  x and  -C  is  given  by 
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(A-24) 


f (x, 


x)  dxd-e. 


Substituting  Eqs.  (A-12)  and  (A-23)  into  Eq.  (A-24), 


dxd£  , (A-25) 


Integrating  with  respect  to  I gives 

0 ' 


w \ 

V 


dx 


(A-26) 


f(x,  i)\ 


2^ 

LW 


f '■  2*^x,  f, 

*'i=0  *'x=0  \ 


0) 

cos  — 

u 


As  in  the  previous  section,  let  us  say  that  the  antenna  is  removed  to  a 
distance  large  compared  to  W from  the  charged  dielectric  surface  on 
which  streamers  are  occurring,  and  that  t is  constant  over  the  dielectric. 
With  this  restriction,  Eq.  (A-26)  becomes 


f(x,  f) 


2 


which,  upon  integration  and  simplification,  becomes 


If (x,  i)  1^  = 2’Jf^ 


1 - 


W -[ 
sin  — ii) 
2v 

W 

2v 


(A-27) 


(A-28) 


3.  Streamer-Noise-Source  Spectra 

The  results  of  the  statistical  analyses  of  streamers  are  given  in 
Table  A-1.  A simple  presentation  of  formulas  in  a table  provides  little 
insight  into  the  way  sti’eamer  noise  characteristics  are  affected  by  the 
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significant  physical  changes  considered  in  the  analysis.  Accordingly, 

Eq.  (A-3)  was  evaluated  for  the  three  geometries  considered  in  Figure  A-1 
and  Table  A-1  for  charging  conditions  such  that  the  total  dc  current 

Table  A-1 

RESULTS  OF  STREAMER  STATISTICAL  CALCULATIONS 


discharged  by  streamers  is  I = 1 uA,  which  corresponds  to  a streamer 
-1 

PRF  V = 666  s . The  values,  chosen  for  the  other  parameters  are  as 
follows : 
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i = 1 m 

-2 

m = 10.6  m 
D = 0.15  m 
Y = tan  0 = 0.58 
W = 0.15  m. 

Since  t = 1 and  I = 1 uA,  the  data  are  In  effect  normalized  to  univy 
coupling  and  unity  streamer  discharge  current.  (To  apply  the  data  to 
other  charging  and  coupling  conditions  it  should  be  borne  in  mind  that 
the  noise  spectral  density  varies  directly  with  the  coupling  factor  ilr  and 
as  •)  The  selected  parameter  values  correspond  to  the  range  of  the 

physical  dimensions  and  couplings  for  the  nose  radcsne  used  in  the  B-47 
flight  tests  of  Ref.  7,  as  illustrated  in  Figure  A-l(a).  The  resulting 
no..‘e  spectra  are  shown  in  Figure  A-5.  It  is  evident  that  in  spite  of  the 
great  disparity  in  the  formulations  of  |f(x,  X)|^  in  Table  A-1  the  three 
spectra  are  almost  identical  in  the  frequency  range  considered.  This 
result  is  significant  in  that  it  means  that  details  of  the  geometry  of  the 
dielectric  region  and  details  of  antenna  location  do  not  noticeably  affect 
the  form  or  magnitude  of  the  resulting  noise  spectrum.  Thus,  we  can  be 
quite  loose  in  defining  the  precise  streamer  geometry  and  still  obtain 
valid  noise  data. 

To  further  explore  the  effects  on  the  noise  spectrum  of  varying 

streamer  geometry,  calculations  were  carried  out  to  determine  the  noise 

spectra  generated  by  a given  charging  current  (I  = 1 .lA , v = 666  s ) 

dc 

on  regions  of  a given  shape  but  of  varying  size.  The  results  of  such 
calculations  for  a triangular  region  are  shown  in  Figure  A-6  and  for  a 
rectangular  region  in  Figure  A-7.  It  is  evident  from  these  two  figures 
that  increasing  tne  size  of  the  region  (this  increases  streamer  length) 
has  the  effect  of  increasing  the  low-frequency  content  of  the  noise 
spectrum.  This  is  in  keeping  with  the  results  shown  in  Figure  4 where 
it  is  Indicated  that  increasing  the  length  of  the  streamers  has  the  effect 
of  stretching  the  induced-current  pulse,  thereby  increasing  its  low- 
frequency  content. 
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NOISE  SPECTRAL  DENSITY.  V G (UJ)  — amperes/  ^radian 


CONFIGURATION 
FIGURE  A-lla) 


CONFIGURATION 
FIGURE  A-Kc)  / 
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FIGURE  A-l(b) 
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FIGURE  A-5  STREAMER-NOISE-SOURCE  SPECTRA  FOR  THE  PHYSICAL  CONFIGURATION 
OF  FIGURE  A-1 
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Appendix  B 

METHOD-OF-MOMENTS  CALCULATION  OF  COUPLING  TO  KC-135  TAIL-CAP 


77 


Appendix  B 


METHCD-OF-MOMENTS  CALCULATION  OF  COUPLING  TO  KC-135  TAIL-CAP 
1 . General 

In  connection  with  SRI's  general  internal  development  and  upgrading 
of  techniques  for  the  solution  of  electromagnetic  problems,  two  method- 
of-moments  computer  codes  were  being  investigated  concurrently  with  the 
program  reported  here,  for  use  on  the  SRI  CDC  6400  digital  computer. 

One  of  these  codes  was  debugged  and  successfully  running  shortly  before 
the  end  of  the  present  research  program.  This  code  was  capable  of 
solving  for  the  currents  and  fields  in  the  vicinity  of  a structure  made 
up  of  a limited  number  of  cylindrical  elements.  Discussions  with  the 
personnel  involved  with  the  internal  SRI  research  program  indicated  that 
the  computer  code  they  had  prepared  could  be  used  very  simply  to  carry 
out  a preliminary  calculation  to  investigate  the  feasibility  of  employ- 
ing computer  methods  to  determine  the  electromagnetic  coupling  between 
noise  sources  and  antennas  on  aircraft.  Accordingly,  provisions  were 
made  on  the  present  program  to  approximate  the  KC-135  aircraft  by  a 
series  of  cylinders  and  to  carry  out  the  calculation  of  the  coupling 
between  the  tail-cap  antenna  and  the  wing-tip  and  elevator-tip  regions. 

It  was  recognized  that  the  approximation  would  be  very  crude,  and 
that  precise  agreement  of  numerical  magnitudes  of  the  computed  couplings 
with  measured  data  would  be  fortuitous.  It  was  felt,  however,  that  the 
degree  of  a^icement  between  measured  and  computed  resonances  would  be 
indication  of  the  ultimate  Wi>ri-,abil ity  of  the  analytical  method. 


Preceding  page  blank 
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2.  Coupling  Calculation 


For  the  calculation,  the  KC-135  aircraft  was  approximated  by  the 
structure  of  0. 5-m-diameter  cylinders  shown  in  Figure  B-1.  It  is 


evident  that  such  a mathematical  model  approximates  only  the  grossest 
features  of  the  airci’aft.  For  example,  the  wings  and  empennage 
structures  attach  to  the  oroper  positions  (jn  the  aircraft  fuselage,  and 
are  oriented  at  the  correct  angle  to  the  fuselage,  but  the  cylinder 
diameter  is  generally  nor  properly  scaled.  In  particular,  no  effort  is 
made  to  scale  the  taper  of  the  airfoils.  The  modeling  was  deliberately 
made  this  simple  in  order  to  ret  ,rict  computer  reeuirements  to  manageable 
levels.  The  restrictions  on  the  modeling,  of  course,  irean  that  although 
the  gross  currents  on  the  aircralt  may  be  accurately  modeled,  (.ne  should 
not  expect  accurate  estimates  of  currents  and  fields  near  the  tips  of 
the  airfoil  extremities. 


80 


The  way  in  which  the  tail-cap  antenna  was  modeled  is  also  evident 
in  Figure  B-1.  The  top  half  of  the  vertical  fin  was  considered  to  be 
isolated  from  the  aircraft  and  driven  by  a 1-volt  generator.  Obviously 
the  size  of  such  an  antenna  is  much  larger  than  that  actually  used  on 
the  KC-135  aircraft,  so  we  should  expect  the  calculated  coupling  fields 
to  be  higher  than  the  measured  ones. 

Calculations  were  made  of  the  currents  flowing  on  the  various 
members  and  the  fields  in  their  vicinity  as  the  result  of  the  1-vol c 
excitation  of  the  antenna.  Of  most  interest  here  is  the  magnitude  of 
the  electric  field  in  the  vicinity  of  the  airfoil  tips.  Calculations 
were  made  of  the  field  intensity  0.5  m inboard  from  the  "airfoil”  tip, 
and  0.01  r above  the  surface.  The  magnitude  of  this  normalized  field 
is  by  definition  the  magnitude  of  the  coupling  to  this  point,  t = E/V, 
as  defined  in  Section  II-C  of  the  main  text.  The  results  of  the  airfoil- 
tip  field  calculations  are  shown  in  Figure  B-2  along  with  experimentally 
measured  data  points  from  the  corona  coupling  studies  of  Ref.  4.  It 
should  be  noted  that  for  both  the  wing-tip  and  the  elevator  tip,  the 
magnitudes  of  the  coupling  are  in  remarkably  good  agreement  if  one 
considers  the  fact  that  no  effort  was  made  to  accurately  model  the 
details  of  either  the  antenna  or  the  airfoil  tips. 

To  better  show  the  degree  to  which  the  computation  accurately 
predicted  aircraft  resonances,  the  magnitudes  of  the  calculated  couplings 
were  adjusted  in  Figure  B-3  lor  best  tit  to  the  experimental  data.  In 
Figure  B-3(a) , the  calculated  coupling  curve  duplicates  the  measured 
resonances  with  remarkable  precision.  Particularly  significant  is  the 
coupling  minimum  at  f - 2 MHz  predicted  by  the  calculations.  When  the 
experimental  data  were  first  obtained  in  the  work  of  Ref.  4,  it  was  not 
clear  whether  the  coupling  minimum  at  f = 2 MHz  should  be  considered  as 
being  real,  or  whether  it  should  be  construed  as  a bad  data  point.  When 
the  data  were  presented  (see  Figure  10  of  the  present  report),  the 
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(b)  COUPLING  TO  TAIL  CAP  FROM  WING  TIP 
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FIGURE  B-2  COMPUTED  AND  EXPERIMENTAL  COUPLING  DATA 
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(•)  COUPLING  TO  TAIL  CAP  FROM  ELEVATOR  TIP 
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FIGURE  B-3  ADJUSTED  CALCULATIONS  COMPARED  TO  MEASURED  COUPLING 
DATA 
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minimum  was  treated  as  real.  It  is  heartening  that  the  rough  computer 
calculations  also  indicate  that  the  coupling  minimum  is  real. 

In  Figure  B-3(b)  the  agreement  is  good  between  calculated  and 
experimental  data  at  low  frequencies.  Above  about  f = 3MHz,  this 
agreement  gets  progressively  worse.  The  magnitudes  of  the  coupling 
variations  do  not  agree  as  well , and  the  frequencies  of  the  maxima  and 
minima  disagree.  Why  the  agreement  is  not  better  is  not  clear  at  this 
time.  Intuitively,  one  would  expect  that  the  majoi  resonances  would  be 
reproduced  even  with  a model  as  crude  as  the  one  illustrated  in 
Figure  B-1 . Unfortunately,  time  and  funds  did  not  permit  experimentation 
with  variation  of  parameters  such  as  cylinder  radius  in  the  calculations 
in  order  to  investigate  the  degree  to  which  the  accuracy  of  computer 
output  is  affected  by  the  simplifications  made  in  the  mathematical 
modeling  of  the  vehicle.  Such  checks  and  experiments  are  essential  to 
final  judgments  regarding  the  viability  of  the  computer  method. 

3 . Brief  Description  of  Method-of-Moments  Computer  Code  Used 

The  basic  computer  code  used  for  the  method-of-moments  calculations 
was  patterned  after  that  developed  by  Syracuse  University  for  AFCRL  and 
described  in  Ref.  19.  To  use  this  code,  one  must  divide  the  structure 
under  investigation  into  a series  of  cylinders  that  are  constrained  to 
run  parallel  to  one  another  at  the  intersections  as  shown  in  Figure  B-1. 
Each  individual  cylindrical  conductor,  in  turn,  is  divided  into  a series 
of  shorter  cylindrical  elements.  01  course,  the  larger  the  number  of 
cylinders  and  subelements  i>ne  uses  in  describing  the  aircraft,  the  more 
accurate  the  simulation.  Unfortunately,  however,  the  computer  core  and 
time  required  to  carry  out  the  calculation  varies  roughly  as  the  s<iuare 
of  the  number  of  elements  employed.  Tims,  it  is  necessary  to  restrict 
the  subdivision  to  the  minimum  number  of  elements  adequate  to  produce 
the  required  detail  in  the  output.  Obviously,  the  degree  of  subdivision 
necessary  increases  as  the  highest  frequency  of  interest  is  increased. 
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FIGURE  B-4  TECHNIQUE  FOR  TREATING  JOINTS  IN 
METHOD-OF-MOMENTS  CALCULATION 

In  the  prelim inai'y  calculations  periormed  on  the  present  program, 

the  aircraft  was  broken  down  into  a total  of  89  pieces.  This  made  it 

possible  to  carry  out  the  minimal  calcuJatlons  ol  the  present  program 

Oh  the  SRI  CDC  6400  digital  computer  with  42,300^  bits  ol  the  core 

memory  of  150,000  bits.  Once  the  program  was  debugged,  it  was  found 
H 

that  this  run  requirecl  61.6  s of  computer  time. 

From  the  foregoing  paragraph,  it  is  evident  that  the  method-of- 
moinents  calculation  extended  to  more  tietailed  configurations  can  tie  up 
a large  traction  of  the  total  core  of  a computer.  Accordingly,  it  is 
not  clear  that  sucli  a calculation  should  be  made  part  ol  the  Air  Fore’ 
Intrasystem  EMC  Analysis  Program.  It  may  be  prclerable  to  generate  tlie 
coupling  data  separately  for  a number  of  typical  configurations  and 
simply  read  the  data  into  PSTAT  as  was  done  here  in  connection  with  the 
experimentally  generated  coupling  data  of  Figure  6. 


85 


ACKNOWLEDGMENT 


The  authors  are  grateful  to  Mr.  S.  Dairiki,  also  of  SRI,  who 
carried  out  the  method-of-moments  calculations  discussed  in  Appendix  B 


Preceding  page  blank 


87 


REFERENCES 


H.  G.  Hucke,  "Precipitation  Static  Ini'-  lerence,  " Proc.  IRE,  Vol . 27, 
No.  5 (May  19:i9) . 

R.  C.  Ayers  ana  J.  0.  Jarrard,  "Aircraft  Precipitation  Static 
Investigation,"  "ontract  W 33-106  SC-70,  Trans-World  Airlines,  Inc. 
(August  1944). 

R.  Gunn  et  al . , "Army-Navy  Precipitation  Sialic  Project,"  Proc.  IRE, 
Vol.  34,  Nos.  4 and  5 (1946). 

R.  L.  Tanner  and  J.  E.  Nanevici;,  "Precipitation  Charging  and  Corona- 
Generated  Inter fere.ace  in  Aircraft,"  AFCRL  336,  Tech.  Report  73, 
Contract  AF  19(604)-3458 , SRi  Project  2494,  Stanford  Research 
Institute,  Menlo  Park,  Calif.  (April  1961),  AD-261  029. 

R.  L.  Tanner  and  J.  E.  Nancvicz,  "An  Analysis  of  Corona-Generated 
Interferenr.'  in  Aircraft,"  Proc.  IEEE,  Vol.  52,  No.  1,  pp.  44-52 
( Januai y 1964) . 

R.  L.  Tanner  and  J.  E.  N&nevirz,  "Radio  Noise  Generated  on  Aircraft 
Surfaces,"  Final  Report,  Contract  AF  33(616)-2761 , SRI  Project  1267, 
Stanford  Research  Institute,  Menlo  Parn.,  Calif.  (September  1956). 

J.  E.  Nanevicz,  "A  Study  of  Pi’ecipitation-Static  Noise  Gi‘ni>ration 
In  Aircraft  Canopy  Antennas,"  Tecli.  Report  62,  Contract  AF  19(604)- 
1296,  SRI  Project  1197,  Stanford  Research  Institute,  Menlo  Park, 
Calif.  (December  1957). 

R.  L.  Tanner,  "Radio  Interference  from  Corona  Discharges,"  Tech. 
Report  37,  Contract  AF  19(604)-266,  SRI  Protect  591,  Stanford 
Research  Institute,  Menlo  Park,  Calif.  (April  1933). 

A.  L.  Hiebert  and  S.  A.  Scharff,  "An  Electromagnetic  Compatibility 
Program  for  the  1970's — Part  I Intrasystem  KMC  Simulation,  Analysis, 
and  Control,"  WN-7256-1 -PR  , RAND  Corporation,  S.inta  Monica,  Cal  it. 
(August  1B71 ) . 


Preceding  page  blank 


89 


10.  U.S.  Standard  Atmosphere,  1962,  National  Aeronautics  a.id  Space 
Administration,  Washington,  D.C.  20402. 

11.  J.  E.  Nanevicz,  E.  F.  Vance,  W.  C.  Wadsworth,  and  J.  A.  Martin, 

"Low  Altitude  Long-Range  All-Weather  Vehicle  Interference  Investi- 
gation; Part  I:  Laboratory  Studies,”  AFAL-TR-65-239 , Part  I, 

Contract  AF  33(615)-1934 , SRI  Project  5082,  Stani  >rd  Research 
Institute,  Menlo  Park,  Calif,  (August  1965). 

12.  A.  Vassiliadis,  "a  Study  of  Corona  Discharge  Noise  in  Aircraft 
Antennas,"  Tech.  Report  70,  Contract  AF  19(604)-3458,  SRI  Project 
2494,  Stanford  Research  Institute,  Menlo  Park,  Calif.  (July  1960). 

13.  J.  T.  Bolljahn  and  R.  F.  Reese,  "Electrically  Small  Antennas  and  the 
Low-Frequency  Aircraft  Antenna  Problem,"  IRE  Trans. , Vol.  PGAP-1 , 

No.  2 (October  1953). 

14.  J.  E.  Nanevicz,  Flight  Evaluation  of  Induced-Noise  Mechanisms  i>n 
High-Speed  Aircraft  ' AFAL-TR-73-317 , Final  Report,  Contract  F33615- 
68-C-1359,  SRI  Project  6091,  Stanford  Research  Institute,  Menlo  Park, 
Calif,  (October  1973), 

15.  R.  J.  Brun  and  R,  G.  Dorsch,  "impingement  ol  Water  Droplets  on  an 
Ellipsoid  of  Fineness  Ratio  10  in  Axis\Tnmetric  Flow,"  NACA  Technical 
Note  3147,  NACA,  Washington,  D.C.  Oiay  195>). 

16.  R.  G.  D<>'sch,  R,  J.  Brun  and  J.  L,  Gregg,  "impingement  of  Water 
Droplets  >n  an  Ellipsoid  of  Fineness  Ratio  5 in  A.xisymmetric  Flow" 

NACA  Techi ical  Note  3099,  NACA.  Washington,  D.C.  (March  1954). 

17.  R.  J.  Brun,  H.  M.  Gallegher,  and  1).  C.  Vogt,  "impingement  o(  Water 
Droplets  on  NACA  65-208  and  65-212  Airloils  at  1 Angle  o(  Attack," 
NACA  Technical  .Note  2932,  Washington,  D.C.  (May  1953), 

18.  A,  C.  Giiibei’t,  E.  J.tnssen,  and  W.  M,  Robbins,  "Detei’mina  t ion  of 
Rate,  Area,  and  Distribution  o|  Water  Di'o'ps  on  Various  Airloils  I i-om 
Trajectories  OI)tain<‘d  on  tlie  D 1 1 ! eren  t ia  1 .Analvzi-i'  ' NACA  RM  No.  9A05  . 
I'niversity  of  California,  Berkeley,  C.i  1 i 1 . (Fid}rua*v  16,  1919). 

19.  D.  C.  Kuo,  B.  .1.  Strait,  ' lmpri>ved  Drogi'ams  loi-  Analv'^is  ot  Radiation 
and  ScattiM'ing  by  Configurations  ol  .Art'i  v ; ,ir  1 1 y Bent  Thin  Wii'es, 
AFCRL-72-0051  , Contiact  No.  r 1 962K-6«-C-IJl  80 , Proiecl  No,  5635, 
Scientllic  Ri-poi’t  No,  15,  Dtpartment  >t  Electrical  and  Computer 
Engineering.  Syracuse  I'niversity,  Syraiuse,  N.Y,  (.lanuary  1972), 


90 


